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CHAPTER I 



INTRODUCTION 

In [21] Payne and Weinbeger presented a method for 
determining pointwise bounds for the solution of Dirichlet 
and mixed type boundary value problems for arbitrary second 
order elliptic equations and rather general domains. In 
the present paper the w,riter will obtain similar results in 
the parabolic case for Dirichlet and mixed type boundary 
value problems where the given data is a function of the 
domain for some problems while for others the data may depend 
on both the solution and the domain of the problem. 

The methods of this paper yield error estimates in 
terms of the integral of the square of known functions and 
thus the results lend themselves to improvement by means 
of the Rayleigh-Ritz technique when the given data is a 
function of domain (the linear case). In such cases it is 
an unfortunate fact that, even if the set of functions to 
which one applies the Rayleigh-Ritz technique is complete 
in the function space of the solution, the actual solution 
of the problem can not be obtained using our present-day 
computers due to the machine creation of large round off 
errors. It is quite possible however that future 
generations of these computing devices will eliminate 
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such errors and thus make the techniques of this paper of 
greater practical importance. In this regard we note that 



2 



the current electronic computing machines ha.ve made of 
practical interest problem solving techniques which were 
previously only of theoretical interest since feasible uses 
were not apparent at the time of their evolution. 

For many of the problems considered in this thesis, 
the existence of a solution has not been established under 
the general conditions necessary for the existence of the 
pointwise bound. When certain addition requirements or 
modifications are placed on the boundary value problem, by 
referring to the literature [l] , |j?] » jj?] to [l 5 ] ? [ 22 ^j and 

[23] we note that the existence of the solution is guaran- 
teed. 

As was pointed out above, the estimates of this paper 
are in terms of the integrals of known functions. For 
another method of estimating the solution function of 
boundary value problems, the reader is referred to the 
excellent works of Fichera [ 6 j to |j 3 ^ where the solution 
function is bounded in terms of the maximum norm. Using 
the results of Fichera as a starting point, Lieberstein [19] 
derived a numerical procedure by which error bounds were 
obtained using a digital computer. 

In Chapter II we obtain some preliminary results 
which shall form the basis of many of the conclusions of 
the following chapters. The identity of Hormander [l 6 ] , 
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which was used by Payne and Weinberger ^2l] „ is modified in 
order to allow its use in a Dirichlet type parabolic problem 
and functions analogous to (4.1) and (4»l8) of [2l] are 
introduced and investigated. Chapter III concerns itself 
with non-linear normally parabolic problems of both the 
Dirichlet and mixed types. For the purposes of this thesis 
a normally parabolic differential operator has the form 




where the matrix is positive definite and Z is strictly 

positive. We shall see that the desired estimates are 
obtained by multiplying the differential operator by a 
parametrix function and then using Green’s theorem to 
generate the solution at the point in question. The unknown 
terms in the resulting expression are estimated as a con- 
sequence of the results of Chapter II and by further use of 
the divergence theorem. In Chapter IV the problems considered 
are degenerate in character, thus the matrix a^J becomes 
positive semi-definite for some problems while the function 
Z need only be non-negative for others. The next to the last 
chapter of the paper outlines a method by which solution 
bounds may be generated for a class of elliptic problems 
which were not considered by Payne and Weinberger in their 
investigation of such problems [23] . In Chapter VI the 
results of the thesis are applied to obtain bounds for the 
solution of a specific boundary value problem. 



CHAPTER II 



PRELIMINARY RESULTS 



1. General Definitions 

Let V be an open (N + 1) dimensional domain of real 
variables (x, 7} = (x\ x^ . . . x N , y) which is bounded by 
two hyperplanes, y = o and y = y o >o, and a surface S lying 
between these planes. As usual we denote the closure of V 
by V and assume that the divergence theorem holds in V. The 
reader is referred to Chapter IV of £ 18 ] for information 
relating to the most general regions in which this theorem 
is valid. Let 



2.1 



Dty) = {(x,y) y = * 7 , (x,y) 6 v - sj 



then 



2.2 V = vUD(o)UD(y 0 )US = VUD(o)UD(y 0 )U S . 



On the boundary of V the Euclidean outward normal is denoted 

by 



2.3 



(nj_,n ) = (r^, n 2 . . . n N , n ) 



where 

2.4 



O N 9 
n y + 2>1 = 1 • 
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Since V is a region in which the divergence theorem is valid, 
it is clear that 2.3 is (uniquely) defined almost everywhere 
on the boundary of V. In addition, unless otherwise indicated, 
we require that S be such that 

N 

X 2 

n, - m 1 >o 

S i 1 

for all (x,y)6S at which the normal vector is defined. 

A function Z - Z(xy)is said to be piecewise continuous 
in V if V may be divided into a finite number of subregions 
in each of which the function is continuous and has a finite 
limit as the boundary is approached. The above is the usual 
definition of piecewise continuity; the following non- 
standard definitions shall also be employed since they 
greatly simplify the presentation of the problem. When a 
function Z(x,y) is said to be continuously (piecewise con- 
tinuously) differentiable In x for (x,y)6V then Z is con- 
tinuous on D(y) for o<y<y 0 and is piecewise continuous in V 

\ rr 

while the i = 1, 2 . . . N, are continuous (piecewise 

continuous) on D(y) for o<y<y Q and piecewise continuous in 

V. Thus if Z is twice continuously differentiable in x and 

piecewise continuously differentiable In y for (x,y)6V then 

Z is continuous in V, the , i = 1, 2 . . . N are con- 

^x 1 

tinuous on D(y) for o<y<y Q and piecewise continuous in V, 
the ^1. • , ii = 1, 2 . . . N, are continuous on D(y) 

^ a 

for o<y<y rt and piecewise continuous in V, and finally 

° OJ 



I 
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is piecewise continuous over the entire volume. 

2. The Parametrix of the Parabolic Problem 

For (x,y)6V let the following differential operator be 

defined 



where the symbol indicates partial differentiation with 



are piecewise continuously differentiable in x for (x,y)6 V 
and are such that 



for any numbers (§•]_, • . • f N ) and all (x,y)6V with equality 
on the left holding if and only if = o for i=l, 2 . . . N. 
The function Z = Z(x,y) is piecewise continuously differentiable 
in y for (x,y)€V. For the point p = (x 0 , y 0 )£ D(y 0 ) 5 at which 
a bound for the solution of certain boundary value problems 
shall later be obtained, we assume the functions a*J and Z 
satisfy a Lipschitz condition for (x,y)6V; that is, we 
suppose there exist positive numbers and A such that 



2.6 




respect to x*, the summation convention is used throughout 
and the components of the symmetric matrix a^*J = a i J(x,y) 



2.7 




2.9 



2.8 
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for (x,y)ev. 



The parametrix function Op is defined with respect 



to the point p by 



2.10 # = 

P L 



for 



4TT (y^-y) 



U 
2 



N-2 

Z(P) 2 VaTp) exp |j-Z(p)p 2 ^{y o -.y))‘ 



2.11 p - (x 0 , y 0 )£D(y 0 ), (x,y)6V, and Z(p)>o 



where 



2.12 /O 2 ~ a^(p)(x r - x^) (x s - x®) 

and a(p) is the determinant of the inverse matrix a^j(p). 
On the surface D(y ( ) - |pj , the parametrix and its deriva- 
tives are defined as the limits of these functions as 
(x^y) -» (x,y 0 ) with (x ? y )6V and x^x 0 . 

Theorem 2.1 . The parametrix has the following 
properties z 

a) At all points (x 5 y)6V - , the parametrix is 

twice continuously differentiable <> 

b) The function 

2.13 (y Q -y) 2 exp (y D -y) 1 z(p)p 2 J(V p ) 2 



is integrable over V for S and A. satisfying 






'$<z 



and 



o . 
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c) For* any bounded function ^ which is continuous 
on Vl/D(y rt ) the parametrix is such that 



2-15 # 5 *. = nt)= 

y<% * 

Proofs From the definition of it is obvious that part 
(a) is satisfied. To establish part (b) of the theorem we 
note that for c^,i = 2, 3 • • • appropriately defined (the 
constant c^ is reserved for later use) 

2 - 16 ,j (tpf Jlv 



i c a fJS ('tr-%) i+ ’ 1 { + (* )f\‘* (’i i 'j //'" ' 



+ 



( 1 ^) + ^ * (f) * V 



r* # r r £ f **/ 

L L J te* c*-*> c 



^3 + 



V& W/ 0 



„ *-/, -tf + A 

+ V 1 ' 




9 






. „ A r r £ r -if -2 +2 v-/ -tf+A 

^ C 3 + M {, J / K a + C 7^ 0*<f£>* 



Z/+.S - V 4* A ^^5 «. tt-?4*sLl 

+%* <* T >* + ** ( ft .-*') 1 }’ 



*^['(i~$) 4 ($^6 ('&>■'$) ] J . r <. dojA ^ 

1 

where 

2 . 1 ? r " “ 51 ( xl ” xl ) 

~~ o 

1 

2.18 g 1 i ^(x k -x k )+s 1 ^(y r ,"y)“ aij °( x »7)“ aij *(p) i s j=l s 2 o . o N 

*c o y o 



with 



2.19 




2.20 

with 



h 1 (x i -xj) + h y (y Q -y) = S(x,y) - 2(p) 



i-l 



a 



• t * t / 



bJ 



2.21 
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Now since the a 1 ^ and Z satisfy 2.8 and 2.9, the functions 

gfj. h. and h are bounded over V and thus the constants 

& k b y i y 

on the right of 2„l6 exist. By integration by parts with 
respect to r and then, y, it is easy to establish (using 2.ll|.) 
that each of the integrals on the right of 2.l6 is bounded 
and thus part (b) of the theorem is proved. 

To establish the validity of 2.15 we suppose that 
£ Q > o is chosen so that for o<C<i Q we have 



2. 22 



l u , If 

"r'j- oft) 


'TZ'XpoU = J 


(va)eY 


Cv^ev 


t -^Lv*v 


j/ . . 




( x ,^) 6 V 


/ 


fjsWO 


r 

J J 







J 2 Xa tis 

0 ' P 






Since 2.22 is true for every t such that o<£<£ 0 and since 
Z and t are continuous at p6D(y Q ), we have that 

iU, if 

D<$ 1 

O',*) 



6-?o ^ 0 

tV$€V 



2.23 



11 



^ /V~ 

£-*0 ^“"^o ^ 0 



60 > o 

N 



fm/sp, ^ v 

{M| fv-<)^ £ 5 } (Ml T^-<fr £*} 



a hj 



'/'(f) ? (f ) ” ( V 77) ® >^v*v ^w- J" S [‘2(f) ^ f(^4 ] 



£-»o ^0 % 



LO 



= nf)^(f)^ l (‘/rr)" ? {^ J C-O£*’iy<«0 ■ ^ + if 

£■>0 y>^o % ' ' x ' J 

*- »/ 

+ -fciOV Iv J" “ft) ^ ^Uj 
r** ^ d 

+ w I 2 . (a/ - 2 ) /i/ (^b'^ 2, -®^[" 2 (p)^(^(#r / ' 3 r)'} 



= ft^tL-f ) 1 ' 1 (^TT ) 1 { i^vv. ii^v. J J &(a/-z) As N '\y-y) 

•^p[- 8 (p)A?{v(- 3 a -.^> j <£„, <Au>} 

= Wp) «(*)*' W' 2 " {|^ £ fe-zXirW )*' 

' ('-j=“ ->) 2 a -®1»)0 [■*(1»)^ 2 (^V’S r 7} dn.Ju}. 



■-**' 



Hence if N is even we obtain 



2-2 ^ iJw fl d-S - 

T% D <^ 
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,L rs-z 



= nf)Z(f)(VTrJZjL^jL^ J fy — )-/. 

£-» o ^ 0 V ' k> ^ 

[" KfO A/ cLn^ cLuJ 

~~^Cf) {W) j if ~(^)/_^ -?( J o)^(^-^ / 

f ^Cf)(V7T)’^ iw {y ¥ (¥)/jia)] 

£ ->° 2^0 J 

-£ r(f) cj^ 



■c/,UP 



where in the last equation it is clear that 6^ is the 
surface of the N“ dimensional unit sphere and P is the 
well known gamma function,, On the other hand 
for N odd we write using 2.23 



2.2^ 



W f^}) = ‘/'(fa) 2 (f) 2 (V77) *ilvu ^ J 

Ofy) ' £-> 0 0 



‘ft * )( 2 ) il A^iu;, 



•To evaluate the integral on the right of the last equation 
it is necessary to make the following change of variables. 
Let 



2.26 © =Z (f)/2%V$y 



( 
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then 

2.27 



and 

2.28 



^ = i(f) 0 



I C^r^) 1 Jjaj 



^fjrS *■ ^"k 



ZC&i 

~ Jl^srrO J -- 
^-+^0 0 

t r 00 - ' 

- £ Cf) * J Q * Mjo (- &) de 



($>) e (- e ) de 



* 2c f )^ r ( i ). 



The value of the above integral is obtained by direct 
calculation. Substitution of 2.28 into 2.25 now yields 

2 • 29 ^ S V S = x «f?) w ‘ (¥)(¥) . - i P(- {) o„ 

= { nf >)( n )~* r ( i ) Q „ 

where use has been made of the well-known recurrence 
relation for the gamma function. The desired result is 
now obtained by noting from page 223 of [2 J 

2-30 4 = 27r*[r(z)]’' 



i4 



and substitution of this expression into 2 , 2 lj. and 2.29 yields 



2,31 



JT Js 

Vty) r 






for N either even or odd. 



3 . An Auxiliary Fnnc t 1 on f or the P arabo lic Problem 

Here our aim is to construct a sufficiently smooth 
function which will be employed in the following chapters 
to obtain solution bounds for the parabolic problem. The 
desired function * 7 ^"" will be defined so that 

2.32 -'F < 00 for (x,y ) 6 V 



and 



2.33 J~ (% ) ^ 0 for (x,y) & V 



where J satisfies the conditions of the previous section. 

In addition to the above, it is necessary that have a 
singularity at peD(y ), which is of a slightly lower order 

Q 

than the singularity of the parametrix ^ p at p£D(y Q ). 

The desired function Is constructed by first considering for 
N *1 

2.34 - ( y<r $ f ( ^-0 ] 



and computing for (x,y ) € V 

2.35 *03.-9)'^ 



1 

*w 



a 

a 

h- 



£#.)( VXg.--a) 
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-fal UN-/) l (f>) («) + ^^4?. 

+ (^)2(^)-^(aA/-/)z-2-C^ 4 ^Cv^f ' 

' Hf [- U W - 0 £ C f>) (<5/V(^ 0 - 

- 4 0-')2 (f)^ J (>k k ) a/f) ^„- *) - 4»v-<) 2d»)3^(f)(^-# 

^4 c ^ )2 

-(w) ? (2N-/)i(f)p 2 i-(vw) ( y ^)^ p X x 



+(»)' 2 (2»-/fiCp)^^a^^i»)(>>->‘*)a^(f)(><-></)(^ 0 -^) 

-(««)' W<) «Cfi) l.(/-4)f- <.M)\f) 

•jay-jo [-(2 .n-i) 2Cp) / o a (iS/v , ('j,- ^)) j 
with g 1 ^, g 3 -^, h and h defined by 2.19 and 2.21. 

*c y i y 

For (x,y)€ V and sufficiently close to p, the second term 
in the brackets on the right of the last equation dominates 
the first and third through seventh terms while the eighth 
dominates the remaining terms. Due to the boundedness of 
the g * s and h's and the piecewise continuity of in V, 

it follows that there exist positive numbers g ^ and such 



that 



2.36 



JCf ) ^0 






jcf) <o 



3 °* f *&-3* s 0 «?} 



Let be such that 



2.37 






% W- (*59^* - * <0 



for all (x,y) G {x,y | (x,y}£ V, o $ y Q -y < S^} . We set 



2.38 



^ = «iin { , S 02 j 



and then define the desired auxiliary function as follows in 
terms of non-negative constants 0 C and >6* 






:.39 ^=Sl 



•% (n>- 1 * 3^) 



(v 3 £ V ^ 0 ~ 



+ oc ^ {- (* w - 0 *<f ) <a l [* w ^.-t>]’' 1) 
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The constants A. and 9r appearing in the last equation shall 

be explicitly chosen in the proof of the theorem which 

follows this paragraph. On the surface D(y 0 )-{p} s the 

function ^ and its derivatives are defined as the limits 

of these functions as (x,y) approaches (x,y ) with (x,y) £ V 

and x x ; thus we have 
1 o 



2 »40 ^f(x,Y Q ) - fam. S^"(x, y) x jtx 

y*y 0 



2.41 

2.42 

and 

2.1)3 



^(x.y ) = dim C ?\(x,y) 

* TV ^ VT ^ 



y*y, 



<5 7 5r ?i ] o(x»y o ) = &ivrF, ± Ax,y) 



y->y c 



x 5 tx 0 
X ^ x 



o 



— <X,yo) = /im ^ (X,y) 

y*y 0 ~5"y 



^ x o 



Theorem 2.2 . The function w (o(. and appropriately 
fixed) is twice piecewise continuously differentiable in x 
and piecewise continuously differentiable in y for (x,y)€rV 
except at p = (x c *y 0 ) and is such that 

2.44 (x,y) € V 



and 
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2.45 ^>0 , J(^") <o (x,y)£ jx s y | (x,y) £ V,o < € s y Q - 

where the coefficients of J satisfy the requirements of 
section 2 of this chapter and, in addition, the condition 

2.46 min Z - m ? >o 

V 

Proof: The function --obviously satisfies the continuity 

and differentiability requirements of the theorem. 

Our aim shall now be to establish 2.44 an d 2.45 for 
properly chosen constants oC and J/r . Prom 2 . 36 , 2.37 and 
the definition of W 7 it is clear that the desired in- 
equalities are satisfied for all 0 C? o where o <y 0 -y £ S 0 
and ~ 0 ° For an< ^ 0 ~ * 8 q We com P u ^ e by 

means of 2.39 

2-47 T er) = T C*) * oc ^ { §| (frtf 



+ e}- 



<*f f o- >) nf>y\zN (&-$'} 
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For (x,y) €-V and y^-y >£ we have 

2.48 m = ^ { - (wr -< ( 2 ,.,)e^ <p) (,V.)C^f ' 

- (v w)'' n ' * (tw - ,) *Q>) ^ (f ) C^.-g./ ■ ' 

+ 0)' 2 «** (>< 

■% 

^y's(ifj-i) Z(p)f l ( 3 r'&) i " z 

+ ' a«"'d) i ^(P f ' (^ ' J -') 2 y fjd- ! '(? ^ (g . - '/)')'} 

foL t 1| c^-4 - * ^ (g 0 "4 ' 

■ -e*f {-(^--)&C^ 2 < 8A/( 3o-a)y'l- 

Having calculated the above, we are now in a position to 
choose the non-negative numbers o C and -$*. Since 2„37 
is valid in the domain of 2,i|7s, it is clear that >. o 
may be fixed so that 

249 J(^0£o (x,y)6 {x,y| py> Qi ° i y D -y s £ 0 J 
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and 

2 *50 l( c f)<o (*>y)£ •[ x,y | ^0 >^ 0 ,o < £ < y Q -y < £ q } 

By inspecting 2.48 it is clear that $r ^ o may be chosen so 
that • 

2.^1 !&)<o (x,y) £ { x,y | (x,y) 6 V, y -y > § 1 

^ O O J ' 

In order to establish that the auxiliary function 
is non-negative in V and is strictly positive in every 
closed set of V, we need only consider equation 2 . 39 . 

Subject to the foregoing, it is clear that the 
theorem is established; however, as a final remark, we 
note that for certain differential operators we may choose 

2.52 % (x,y)€ £x,y |(x,y) € V,o < y Q -y £ 

since J(^F) satisfies 2. 49 and 2.5>0. When this condition 
is applicable, one would of course choose cC = o in 2 . 39 . 

The conclusions of Theorem 2.2 shall be employed 
in the next chapter to obtain pointwise bounds for the 
solution of non-linear normally parabolic problems. In 
the first section of Chapter IV the function Z- = Z(x,y) 
is not assumed to be strictly positive and 2.46 is replaced 

by 
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2.53 min 3 > o 

V 



and 



2.54 Z (x,y ) >o 



U,y)€ D(y 0 ) 



Under such an assumption it is obvious that Theorem 2.2 
is not valid since we would in general be unable to 
choose 30 that 2.48 defines a strictly negative number. 
Because of this reason we use the fact that a^J is a 
positive definite matrix and introduce a slight different 
auxiliary function. This new function is much more 
difficult to construct since, besides having a more com- 
plicated form, it is necessary to calculate a^(x^,y) 
for oiy-$y 0 . Thus, even though the revised auxiliary 
function could be used in the linear normally parabolic 
case, the function as given by 2.39 would be easier to 
construct since the inversion of the matrix a-*-<3 is only 
necessary at the point p £ D(y Q ) • 

When conditions 2.53 and 2.54 replace 2.46, we 
define as follows: 

2.55 <zf = (^-^) ^ Mf 
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for (x,y) £ V where 

2.56 /° = a iJ (x Q ,y) (x 1 -^^) (x^-x^) 

For our present definition of we notice that 2.35 is 
valid with 

2.57 (0 replacing 
and 

2.58 a rs (x o ,y) replacing O^fp), 

In the revised form of 2.35 it is also necessary to add 
the following term on the right side of the equation 

-i 

2.59 sf »%-£>**) ■ 

By exactly the same arguments as previously used, we deduce 
that there exist positive numbers 2 £ an< ^ suc ^ 



that 
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2.60 J{GF)$o (x,y)e[x,y - 7 0 ~^^ 2 ^0ll 

and 



2.6l 



J(^) <o (x,y)^{x,y y o"^ i25 oJ< 



Since Z(x,y) is strictly positive on the surface 
D(y Q ) it follows from the^ conditions of the first section 
of this chapter that there exists a positive number £ ^ 
such that 

2 . 6 2 §(*-#*(■ * < a 



for all (x,y) f [x,y |(x,y) 6 V,o < y Q -y< 2S q2 ^ 



We set 



2.63 




and notice that 

2.54 J(c£ ) < o U,y)6 jx,y | (x,y)£V,o$y 0 -y£ 2S Q j; 



and 
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2 . 6 ^ 



J(°C)<o (x,y)6{x,y (x,y) 6 V, o<£ i y o ~y£2S Q ^ 



where 

2.66 ^ = oc (fa-ft) * 20^0 {" O- i) *( p)fo*i ($,-•&)) ’} 



Since 2.62 is true, it is evident that 2.64 an< 3 2.65 are 
valid for all dC> o. By comparing <^T and of , we observe 
that we may choose oC - such that 



2 . 6 ? 



) < o 



(x>y)€{x,y 



(x,yU V,o<y 0 -y<2§ 0 j 



and 

2.68 ) <o (x,y)e{x,y |(x,y) €■ V,o<Ci y 0 “y$ 25 o ^ 



Next we introduce the function ^ given by 



• {- (ms 0 y'/3 (zw-i) /S' 1 } 

4*^ 1" o? 
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where for the present /O is an arbitrary non-negative 
number. To fix this arbitrary number we compute 

2.70 ?(£•) = { s; 1 -c^-^y }■ 

' {- (WS,)' 1 /! (xv-/)s,p) Cx,^) C**-*/) 

"Gw8o)/3 (««-') *tys) <£^< 7 .^ C *. 3 3 ) 

t(wSj 2 rfV* (zv-,)* 

• C*-*.*) 

-(.fy-i)* ^i-(M s 0 )> 

•h J ( -bcC ^ * 

For y o -y ^ 2 and (x,y) £V, the last term on the 

right of the above equation is a strictly negative 
number. Next we consider the following limit 
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p->o 0 

'{-(vs 0 )''^ (iw-o 2-Cp) + %^ j 



By the remark following 2.70 and due to the form of 2.71, 
we deduce that there exist positive numbers ^Q 2 an <5 (3 
such that for 

2.72 p = ^i 

we have 

2.73 J(#")<o (x,y)e{x,y | ^ Y°02 ,o ^ 7 o“ S Q } 

We are now in a position to define the desired 
auxiliary function. As was mentioned above, this function 
will be used when bounds for the solution of a degenerate 
problem are sought. The function is explicitly given by 
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o£ 



(x^)4V . 



(° Yox , °Ys^o' S o 

4., 



• 74 f 1 ■ t V-c^yH-c^^cw-!)^ (i)V 



+ 



(‘fu$ t ,<D''(™-i)Hp){3 l fX +>} 



■ ^p{- (ws^y $ (2 a.-/) 2-c r )^ j 



4d +o( 



(° (°oi j 0/S 3 ' $> 



On the surface D(y 0 )--fpj , the function d" and its 
derivatives are defined as the limits of these functions 
as (x,y) approaches (x,y Q ) with (x,y)6V and x^x q . We 
now show that wi th properly chosen the function 
satisfies the conclusions of the following theorem. 

Theorem 2. 3 . The function‘d”, defined by 2.74> 
is twice piecewise continuously differentiable in x and 
piecewise continuously differentiable in y for (x,y)6 V 
except at p = ( x o’^o^ an ^ suc ^ that 



2.75 


J {^F) * 0 


(x,y) < 


5- V 


2.76 


J(<=F) < 0 


(x,y)€fx,y 


(x,y) G V,o<£ « y Q -y 


and 








2.77 


ma X {-T} 


= M 





where M is a finite number. The coefficients of the 

differential operator J satisfy the requirements of 

section 2 and in addition the condition that the partial 

derivative of CL (x ,y) with respect to y is piecewise 

rs o 

continuous on the line (x Q ,y) £ V. The function S is such 
that 

2.78 min «^o 

V 



and 

2.79 Z(x,y)>o (x,y)eD(y Q ). 

Proof: The continuity and differentiability of are 

obvious for (x,y)e{ x,y|(x,y)GV, oSy Q -y£ % ^ . 
lower part of V we calculate 



For the 
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2.80 



and 

2.81 



&l' At '{ST'-Cyu)"]' 

O'Vi.-h 



P > Po l 

0 i pfc-s* 



Hi =cc-(^r}' 

P'fal 

*t"(7^ 0 ) ^ (i^-i)s(f) ^6v,^ 

'^f {" (*A> SoT 1 /? (2W'j)*(|,)^ l z } -f 

= Awu {^1,?. 

A-VU * 

H >/H 

0 4 ^ K 

Our task is next to select the number <</• so that 'fr' 



satisfies 2.75 and 2 . 76 . Prom the remarks associated 
with equations 2.55 to 2.73 it is clear that these in- 
equalities are satisfied on the set 
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2.82 






For jO >p$2 and °< 7 ^.y 0 " S 0 we compute 

2 . S 3 7fr> { { Osj' 1 




'**# {- Os J’fl O-j)^),^} + + 

+ fe"- f^-^r'] 

-t(vusj)' 1 (iw-i) -eifififil + i } - 

Sj'fi, (zw-i)*(f^3 + 7 (&**). 



In order to determine an appropriate choice of the 
constant JL we note that 



2.8^ 



a i Ja ir (x 0 ,y) (x r -x£)a sj . (x Q ,y) (x s -x^)^ a^ 1 ^) 



2 

02 



for all (x,y) belonging to the domain of definition of 
2 . 83 ; hence it is evident that may be chosen so that 
2.83 defines a strictly negative number. 

With d_ so chosen it is obvious that the finite 
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number M given by 



2.85 



M — max] 
V 




exists since i 7 is bounded everywhere on V except at the 
point p and at this point one has 



2.86 



sdi/YYU ^ o , 

F eV 



In future investigations bearing on the subject 
matter of this thesis, an attempt will be made to replace 

2.79 by 

2.87 E(p) > o . 



Next we shall establish the existence of an 
integral which will often appear in the bounds which are 
developed in the following chapters. 

Theorem 2.1; . The function is such that the 
integral 

2.88 Jv 



exists , 
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Proof: As in the other theorems of this section, the proof 

is completed by discussing two cases. When 2.^6 is satisfied 
it follows from 2.3f? that there exist positive numbers S ^ 
and ^ which are such that 

- 1 . 2 



N+3 



2.89 - 7 ($)>Cy-y) 2 



|z(p) (y 0 -y)+(32N 2 )” 1 (2N-i)z(p) 2 ^> 

' exp { I"- 8 N(y o -y )]" 1 j(2N-l) Z (p) 2 

for ^ and oi y Q -y<S 1 1 & Q . With jO 2 replaced by 

p> 2 it is evident that a similar result follows for Of 

in 



case 2.79 is applicable. Hence we have 

- N+1 

2.90 ^(y 0 -y) ~ 2 ~z(p)expf- [ 8 N(y o -y)J 



-1 r. 



(2N-l)Z(p)^ 

for ^ ^ and o £ & . Next we consider 



2.91 



= lfJ[-jfr)}''[Kn r pv + UflrmJ'inirfjv 



where 



2.92 



v ! ■ {< x >y) P’P -</V oi v 7S ^i}- 



The first integral on the right of 2.91 clearly exists since 
its integrand is bounded over the domain of integration; 
also, since 



(8 n ) ** 1 (2N-1) < ~ 



2.93 



9 
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the second Integral on the right of the aforesaid equation 
is seen to be finite by application of 2«90 and Theorem 2.1 



3* Generalized Form of an Identity of Hormander 

The identity under consideration was derived by 
Hormander £l6j for the purpose of obtaining a uniqueness 
theorem and estimates for normally hyperbolic partial 
differential equations. Essentially the same identity was 
employed by Payne and Weinberger [21J to obtain solution 
bounds for second order elliptic equations. In what follows 
a slightly generalized form of the same identity is derived. 

Let f^(x,y), . . . f^(x,y) be any set of functions 
which are piecewise continuously differentiable in x for 
(x,y) 6 V. By direct differentiation and application of the 
divergence theorem, we have for any function 'f' which is 
continuously differentiable in x for (x,y) 6 V and whose 
second derivative in x is continuous in the interior of a 
finite number of subregions, the sum of which is D(y)» where 
y is such that o < y < y Q (in the above we suppose that D(y) 
is a domain of N dimensional space) . 



2 -* //[fV*- 4* W 



%. </,. rtv.ls. 

* A 
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The functions a 1 3 i,j = 1, 2 . . . N, which 

appear in the last equation, are assumed to be symmetric, 
piecewise continuously differentiable for (x,y)6 V and such 
that 2.7 is satisfied. (For the results of this section the 
piecewise continuity of the derivative of a*^ with respect 
to y is not needed? but, since this condition shall later be 
necessary, it is included here.) 

Let 

2.95 ^ 



and then note that the vector 
2.96 



(T* 0) = (a. J ^ ») 



defined almost everywhere on S is orthogonal to (n^, n^) 
where in 2.96 we have made use of 2. 5 and 2.7* From the 
above it is clear that (*7^» 0 ) is a tangent vector to the 
surfaced, is perpendicular to the vector (0,1) and is 
defined almost everywhere on S • Following closely the 
procedure of section 2 of [21J we renormalize by setting 



2.97 



ctSo) - (rVV 2. o) 



rs 



where a^j is the inverse matrix of a*^. For (x,y)6<5 



we have 



n j n 

where n is defined by this equation. Due to the above, the 
directional derivative on $ , 



2.98 



a i ^t i t*^ = a^n^ 
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2.99 



V£ _ 






" V f K t 



la in a tangential direction which is perpendicular to the 
"y" direction. For (x,y)€ S we notice that 



2.100 



• a 



,*“ 1 v 






V rv — . 



- c i 



v . \t dk 

\ iV -n _ 



C-sk 

= 7 f 



\ ' 



Also 

2.101 



a,// T 7 ,‘> 



/a of 

j; V ^ /M3V/ ^ 









/ ( \ y 



a«A*- 



Hence almost everywhere for (x,y) £ 5 



2.102 



& *j-^v =/n 






kJ*T 



/ ^ 



*/ 

vn. 









Moreover almost everywhere on 5 we have 
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2 - 103 






V" <7 



a..fV* 



nvf 

V " - , 



/w 



j> < . r ^ ^ ^ 

= /rt f - f ^ , 5D 



/.«> ^ , „<i,A^_ A * 

f 5 V + ^ . 



and thus 
2.101+ 

Using 2.102 and 2.104 we rewrite 2.94 as 

S.1M J/-1»\{(|.')‘-«) t )^./< , Hsl]Js 

rrrlr^ ** ^ ^ h a / 

r^“- 'V 'V +f 

1c 

We now choose f , k = 1, 2 . . . N, so that a bound for the 

square of the integral of the conormal type derivative 2.95 

may be obtained from 2.105* Previously we required that 

these auxiliary functions be piecewise continuously differ- 

*? 

entiable in x; we now further assume that f n^ s is bounded 
and has a positive minimum on 5. Such a set of functions 
was considered by Payne and Weinberger in section 3 of |^2lj . 
Prom this reference we note that, if 5 is starshaped in x 
with respect to some line perpendicular to the x plane, then 
we can. choose this line as the x origin and take f^ = x^. 
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Also if V is a domain between two boundaries each of which 
is starshaped in x with respect to the above perpendicular 
line and if V contains a cylindrical shell of radius r 

o 

from this line then we may take f k = x k (r-r ) . For further 

o 

information bearing on this choice of the auxiliary functions, 
the reader is referred to Appendix A of |~17 ]•*'*" With f^ as 
defined above, set 

2.106 m 0 = min n“ 

3 S 




Since a*j is positive definite, there exists a constant 



Ci which satisfies 



'1 
2.107 



[f, k a i3 -f,*a 3k -f, 3 a lk + 



for all (x,y) 6-V. (C.^ is an upper bound for the largest 

eigenvalue 0^ of the coefficient matrix on the left with 

0 o 

respect to a 1 J.) For a crude upper bound for we define 



2.108 



Jr 



ij = 






f 'k a 



f,Ja ik + f “■a 



k n ij 
*k . 



then 



**See also, Bramble, J. and Hubbard, B., Some Higher 
Order Integral Identities with Application to Bounding 
Techniques, in print. 



38 



2.109 






c ' " = 4“ !, f, 



* 









B 



^ a ~ C. 



where 



2 «110 B - max j/~ J ) f 

V C /$ J 



Substitution of 2.106 and 2.107 into 2.105 yields after 
applying the indicated inequality to the resulting expression 



2.111 



S 5 o 

rSU-’f\(W^ *2 SSJf\ («.\\ 3 J-V 



fC. JTJYV,; 



V 

where c£ is an arbitrary positive number. By setting 



2 . 112 



& = 1 



we may write 
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2.113 




.m * 

* v 






where 

2.114 









and clearly 



2.115 




*= 0 



if r = o for (x,y )6 S - In the following chapters great use 
shall be made of 2.113 when the boundary value problem under 
consideration is of the Dirichlet type. 



CHAPTER III 



NON-LINEAR NORMALLY PARABOLIC PROBLEMS 



1. The Non-Linear Mixed Problem 

In this section, bounds for the solution function of 
a ’boundary value problem, which is normally parabolic, will 
be derived where the solution function is initially given 
and thereafter satisfies a non-linear Neumann condition on 
the boundary and a given differential equation in the 
interior. 

We assume that V is an N + 1 dimensional domain with 

boundary D(y )U SUD(o) which is such that the conditions of 
o 

the first section of Chapter II are satisfied. Our goal is 
to determine pointwise bounds for the solution of 

~co< W ) VW<°0 



3.1 



< 



W(x ; o) 



&(*) 






- a *w Jt on,. - ICxj^w) (K$)eS 

^ — 00 < W < OO 

where V/ is continuous in VUD(y Q ), continuously differentiable 
in x for (x,y) 6 V and is such that its second derivative in 
x is continuous in the interior of a finite number of 
subregions, the sum of which is D(^), where y satisfies 
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o <y<y Q . Also we suppose that the first derivative of w in 

y is continuous in the interior of a finite number of 

subintervals s the sum of which is the line (5r s y) for every 

fixed x where (x,y)£V. In addition it is necessary that the 

solution of 3.1 satisfy 
r 

lim w(x,y) » g(x) 

(x>y) ->(x,o) 

(x s y) € V 

3-2 J 

lim a^ J (x,y)w, (x,y) = a 1 J(x»y)w, (x s y) 
(x»y)->(x,y) 1 1 

(x,y)€V A, - .... N , 

(5,y)6S f 

The boundary value problem is assumed to be such that 
the components of the symmetric matrix a*3 = a* 3 (x,y) are 
piecewise continuously differentiable in x and are such 
that 2.7 is satisfied in V. In the closed set V the function 
Z = E(x,y) is piecewise continuously differentiable in y and 
satisfies 

3-3 min Z(x,y) =m > 0 . 

V d 

The boundary data "f , Jl and ^ are assumed to be integrable 
and square integrable over their domains of definition for 
any bounded continuous function w which is continuously 
differentiable in x for (x,y) £ V. In addition f and JL 
are required to satisfy a Lipschitz condition in all but 
their first two arguments. Hence there exists positive 
numbers M , i = 1, 2 . . 



. N such that 
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I f (x, (x, y , w* ,X7 Wg) | • 




3-4 -< 



(^)ev ^ - 00 < v/j 7 w < 00 



* f\ l v/ ( " w zl 



(*, ^) € 6 3 — oO < v/ < oo . 






Except for the fact that H - &(x,y, vy) is required 
to satisfy a Lipschitz condition in w, the problem of this 
section is a generalization of a problem whicn Friedman 
considered in JjLlJ . In that paper the elliptic part of the 
differential equation was the usual Laplace operator, •f 
was identically zero and Ji was a monotone non-increasing 
function of w . Under these restrictions the solution of 
the boundary value problem was shown to exist by means of 
the Schauder Fixed Point Theorem [26]. In a concluding 
remark the author indicated that in a later paper he would 
attempt to extend the existence proof to general second 
order parabolic equations. Whether or not such an extension 
has been carried out is not known. The existence of a 
solution for a similar problem was however established by 
Piskorek [22] and also by Pogorzelski [23]. 

Let p £D(y Q ) be the point at which a bound for the 
solution is desired. We assume that the functions a*3 and 
2 satisfy a Lipschitz condition at p€D(y Q ) for (x,y)£V, 
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that is, we assume that 2.8 and 2. 9 are satisfied. In 
addition it is necessary to suppose that f(x s y,w, vw) is 
bounded in some neighborhood of the point p which is con- 
tained in V for any continuous function W which is continu- 
ously differentiable in x for (x»y)£V. Hence for any such 
function w there exists numbers 8„, /0 t and ^ (w) all 

greater than zero such that 



3-5 



f (x,y,w, v w) $ °Hr{ w) 



for p and y Q - 6 2 <yiy Q . 

The desired bound for the solution is obtained by 
choosing an arbitrary function IP (x,y) which is such that 

^ Ip 

J(iP), Jv and ^ (x,o) approximate f (x,y, ^P ) , Z(x s y, <J>) 

and ^ (x) respectively where we further require that (P be 
twice piecewise continuously differentiable in x and 
piecewise continuously differentiable in y for (x,y)£V. 

Let 



3.6 



^(x,y) = w(x,y) - <p(x, y) 



and compute 



3-7 



r 






(x^)£V 7W \7 4? < oO 

iv ■ i (*, ^,w) -i (*, ^ ip) + L(xj a^) 

(< ; ^) ^63 ~~ c ° ^ ^ 



where 



1 
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3.8 

and 


P(x,y) - f(x,y, IP s V ) - J( p) 


3.9 


n %lp 

L(x,y) = Ji(x,y, p ) - ^3 . 



Due to 3.3 and 3-4 we rna y write 



3.10 

and 


J( V) fMV,. J + F(*,^) 


3.11 


Tv x 1 F| + L (*, ^ ) . 



With TF(y') and ^ defined by 2.6 and 2.10 respectively, 

P 

from the divergence theorem it follows that 



3.12 


XT/{y f jc^)}Jv 



= JJ fyVS - * 

%) PC's) 



where 


JT /) V^> ^ y <Ls + If Uf 

d(g) ? ' D(e) r 


3.13 

and 


V(y) = ‘[(x,/y) (x,a|)6D(^), o</Yj<yJ 


3.14 


S(y) = { (x,/j ) | (x,^ ) 6- 3(nj )-D(^ )> 0 </^<y]. 
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In 3. 

3.15 

which 

3.16 



where 

3-17 



.4 we let y approach y Q and use Theorem 2„1 to obtain 

*<*',*'> = ^ W)} ^ 

+ II{jf r ^ - t 5 ^ } is - ^ 5 + If HpiVeLs 

becomes after substitution of equations 3*7 to 3°H 



lw(x 0 ,^) - %(f>)l 



+ 5f *p '<’)} <4 

\S 

- If f Ty ^ s ~ 2 V ds | 

$ s 0 r 



* i* i^^w***. 

■{ ms f {<fi ( 't n^lsh 



X(f) ty( y O)'$o) + £f Ux^'fy^S’ 

+ JT ^ & ds ~ f dv . 

D(o ) r V(^) / 

The terms on the right of 3«17 are all known; 



however, it Is necessary to prove that the last integral 
converges o For this purpose we write 
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Ifj ^AO F J.V 
7r*% v(^) i 



— iiw. ^ /// x> Fjv f §fj f J.v J 

W ^-v(aa) vfojw^) 

v-v(ftA) v-v^s5 



where 



tiwv I// i* {?(A'j,ie'7ip) -j(tMJv 



3 - 19 V (AA) = {^^>1/°% , 

and on the surface D(y 0 )-{p} it is clear that X is defined 

P 

as indicated in Section 2 of Chapter II. We now consider 

the three integrals on the extreme right of 3.18. The 

first exists since ^ is bounded in the domain of inte- 

P 

gration and f is integrable over V. The second integral 
exists since its integrand is bounded over the domain of 
integration. For the last term we use the fact that 
[f (x,y, <4 , ) -J( *f )] is bounded over due 

3.5 and the definition of J; thus the integral exists since 
"if is easily seen to be positive and integrable over 

-T* 

v(/° 2 , S 2 ). 

From the above remarks it is clear that all the terms 
on the left of 3*l6 except w(x 0 ,y 0 ) are known and hence to 
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obtain a bound for the solution it is only necessary to 
estimate the expressions on the right. To obtain this 
result, we introduce the function T defined by 

3.20 V' = ^exp[-K(y 0 -y)J 

and calculate 



3.21 J(V ) = J( ?exp [-K(y Q -y)]) = J( £ ) exp [-K(y 0 -y)J 
where 

3.22 J(£ ) = J( p)-Af ? 

Let f*, i = 1, 2 . . . N, be the set of auxiliary functions 
which was introduced in Section 4 of Chapter II, that is, 
the f 1 are piecewise continuously differentiable in x for 
(x,y)-6 V and are such that 



3.23 



o < m^ = min f n^. 



With the f* as thus defined we may write 

3.24 I/rUs 



v 






where cC is any positive number and 



3.25 



= max 



{X f 1 ^} 
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3.26 




In arriving at 3.24 we have made use of 



3-27 (fV.f J E E-fV +■ 

’* l#t ‘i 'lr 

which is easily established by expanding 

3.28 E (4‘t '-{* % ,f. 



Prom the divergence theorem it is clear that 



3.29 



UStHVJ? 








Jv + JI5 I %ls 



S 



- Z M $ z zjs z^G^zJis 

V Dfy) D0>) 

+x JT/ i ?^(|f)^v - { 



Using 3«7 and 3° 20 to write 



3.30 m, itl + n l ifj + f<H^)^jb[K(g»-»)] 

and 

3.31 Tv ' + L (*, q) [kty-*)] ^ 



we have by means of the above 



3-32 
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Itfmwdvl* ju\n\muv 

V V 

<», sss 9 + tiU'% v h t ss ^ 

+ 'I 4 XjX '/' olv f ^ fff F Z Jh^p[£K(^-^.)] Jv 

and 

3 ‘ 33 Iff? s$«fsl S + T^SI^S 

s ^ s s 

*■ ST, XT i.^^p.KCg.-'j)] Js 



where cC 

3-34 




and are arbitrary positive numbers and 

-t (M^) 2 . 

/ 



Substitution of 3*32 and 3*33 into 3*29 yields after the 
necessary transposition 



3.35 



0- ftyWv 

s (V f + r«s-)Xf ?Vs + 



where 

3.36 




and 
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3*37 NL = max {-n z}. 

i> s y 

Except for the first and last integrals, the right side of 
3.35 is completely known; the last term may, however, be 
neglected since Z satisfies 3*3* In addition we note that 
the left side of 3*35 has the same form as the right side of 
3.2I4.; and, after a close scrutinization of these inequalities, 
it is evident that they yield a bound for the integral of the 



square of f over S provided the arbitrary positive constants 
oC^ and K are properly chosen. In a particular problem the 
best choice of these arbitrary constants would depend on the 
functions P, G, and L and the M's appearing above. For 
definiteness we assume o and make the following selections. 

(If = o, obvious modifications are made to obtain the same 
result) . Let 



308 

3.39 

3.40 
3.i|l 

and 

3.42 



a 



cC = — 
2 



<*. = 1 



oC = 



a^mi 



°it (M 2 + \ + \ M^) 



4m 



-1 



*3-1 



K = 





+ 



1 

2 



+ (2M 2 -H+M^) 



:(2m 2 4-i+m^)m 1 

2 

a m, 

o 4 
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where we have assumed that 
3-43 M 2 + | + | M^>o. 

If this is not the case, then one need only choose °C ^ so 
that 

'•w cC ), -i 

3.44 M 2 + 2 + 5 M^> 0 

and the same result is now obtained by appropriately defining 
K and °C 2/ Substitution of the above into 3-35 and 3-24 
gives 

3.45 XffSh s (^^+i+n s )'{sj 6 : ‘£^(z^ l! )Js 

where is defined by 3*4^ (assuming 3*43 I 3 valid). We 
use the last inequality and 3.20 to obtain if 

3.46 |j i (2# 2 +i+ n^y Q*i js 

+ SfS F^^zk^Jy 4 X[l 2 <hf (- 2 k, q) tl s] 

and if 0 then 

SSfti i (zfi-H* /$)'{£{<: l ir^*p (tK^yJs 

S No) 9 v y 



3.47 
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The right side of the applicable inequality is known since 
F* L and G are square integrable over their domains of 
definition; hence we may write 

3 . 4 a If v'Sh <13 

3 

where 1$ is a known number depending on the function ^ . 

After having calculated we note that a bound has been 

obtained for the last term on the right of 3.16. 

If a function P p can be determined which is such 
that the conclusions of Theorem 2.1 are valid and in addition 
the condition that 

3.49 J( r p) = 0, 

then it is clear that Pp is a fundamental solution of J and 
the first term on the right of 3«l6 is identically zero. 
Dressel [ 3 ] and [4] established the existence of such a 
solution when the coefficients of the differential operator 
are sufficiently smooth. From these references it is clear 
that even when its existence is guaranteed such a function 
may be quite difficult to construct; hence we shall assume 
that Pp is not known and that the first term on the right 
of 3°l6 must be bounded. To obtain the desired bound we 
introduce the auxiliary function defined by 2.39* Since 
the differential operator of this section is Identical to 
that of Theorem 2.2* we conclude that OC and ir may be 



fixed so that 



3.50 


^ z 0 


(x,j)e v 


3.51 




(x,y) € V 


and 






3.52 


; )<o j ^ > 0 





for every i>0 and y satisfying 0 <££ y 0 -y-y o * Next we 

~ ** 

define the adjoint J of J by writing 

3-53 J( f ) = j( ^ ) - KZ 

Because the singularity of ^ at p is of a slightly 

lower order than the singularity of the parametrix y , it 

P 

follows by the proof of Theorem 2.1 that 

3.51* jL. II 

111’ «3) 

Is^o 

= 0 

since Z is strictly positive and ^ is a non-negative number. 
By means of Green's Theorem we write 

JU U{frW')-V l 5ffi}<LY 

t-»io v 

= 11 ^- 9 l TiUs -SUfFi^ds 

s s 3 

0(o) 



3-5£ 
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where we have used 3-54 and the fact that 

3.56 fS ^ - p 3 - ^ } M 

'i ifj) 

for p £ D(y Q ) • If we assume that K is non-negative, then 
expansion of 3° 55 yields 

3.57 1U HifcWtfC?) - ? l 7(*F) 

v(^). y 

<1. 

+ ^F dv 

S <, 0 [lK^,-$)]cts + x t Sf y^Js 

S 

- If Tv Sftfi fJs 

sS >5 5 d 

fjT^V-a t*f> &Kfy)ds 

where in the right inequality we have for definiteness set 
= 1 and 



I 






I 






3*58 M 6 = max ^ s}. 

By means of 3*7 and 3*21 we have 

3-59 Jfj2 q F9 , ?(P)clv = Z JTT^ p F^fo [Kfy,-*)] <^1/ 

+ £ iff?" ? {f fc^WjVw) -f 6sty *1 V, 

v (l) 

Also 

3-6o X. JT.f '** I Pi • I -fO, y,w, ^w) 

v(*) • 3j 3 

- f (*.•3, * *)l <*)] A 

JIfTP*A 

VCD W 

s (sift, + -Q JflV ?Vv + 0- ffJW i h Jv 

Vty " V (D) * 



where is given by 3*34 and ^ is an arbitrary positive 
number. Let 



3.6l 



M 

OC . = 

5 2a 



then by choosing K so that 



3.62 



M 



K = ?rr~ ) 2M_ + 

^ m p l 1 2a 



+ _l.y = 



K. 



and using 3 * 60 , it is evident that 3»5>7 may be rewritten as 

V% %) 
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With ^ defined on D(y 0 ) - {p} as is indicated in Chapter II 
and with the constants S ^ an d cilosen so that ^ ^ and 

& 2^^o s we ■* aave ^ rom 2«35 for appropriately defined constants 
j 1 — XO $ o • • 13, 



3.64 



hi, { V F [/1 u, -/^)1 -f Z j (*F)} J_\j 

r*ito v-v(-j) ' L 0 J 

- TC+)} J/ 



Jlj U.^h, 
r*. v «A)-%) 



_ A/f l_ 

(V'j) *■ 



_Nj3 

^ e .xho-%) *■ » 



Jt 



tc o (■ v i>' V '] ^ (-[sa/a-^]’ 1 ^-^^])} A 



= 0 



In arriving at the last result we have integrated by parts 
first with respect to ^ and then with respect to y (see the 
proof of part a of Theorem 2.1). Prom the above it follows 
that 

3-65 f[\(V’i)]-P* J&jji'/ $ k . 



Prom Schwarz ' s 



ci t . 

inequality and the fact that V satisfies 3*51 



it follows that 
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3.66 -USS^VF [\( V ^)]^V| 

>-{~m v'jwjv}*- 

-{ fff [-j<&}''(?F) l *r [*K t <v '*)]^} 1 

The factors on the right of 3*66 each exist. The first 
exists since is continuous on V and is hence bounded on V 
while J(^) is easily seen to be integrable over V. For the 

p 

second factor we note that the singularity of at p€ D(y 0 ) 
is less than that of J( ^p) 2 and the result may thus be 
obtained direotly from Theorem 2.1. Due to the above 

3-67 {///-* (<F) 1 

t V F-ty-p 

and then, after substitution of 3*65 and 3*20, we write 

3-68 T(«T) ^j v } k 

< j v f> 



*ty t * sss[-m] 



Because of the form of and 
inequality to obtain 



we use Schwarz ’ s 
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3-69 liU, SSS VT(X*)d\/\ 

Ui) r 

• (i f Jf [-■ W )] [a k 2 (% - $>] J v } * 

+ {% * i") 



where the first factor on the left exists due to Theorem 2.l| 
and the second factor by the remark following inequality 
3.66o It is clear that 3*69 is the desired bound for the 
first term on the right of 3*l6. 

To obtain the desired bound it is now only necessary 
to bound the second term on the right of 3-l6o We write 



3.70 



l^vvv fjjf Xp { -f W ; 7 w)-f (x, ^ V 

<L^ itfxAPi.WI t nVjUv 

4 { /// Jlfft I* M fjftvf 

(7 C/ 0 ^0' 
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In the last Inequality the first factor on the right exists 
since for c^ appropriately defined 

3-71 < c, ¥ t 

Tfo £ f j 

°o N 

for t sufficiently small. The integral in the above 
expression Is seen to exist by means of integration by 
parts as in the proof of part (b) of Theorem 2.1. The 
last factor on the right of 3.70 Is now bounded by recon- 
sidering 3. £7 and 3*59* From these inequalities we obtain 

3.72 (l jjj t&F ?.<,'*» +2KiW'-}Jv 

0 yi, v w) " 1 

A $ < ^ 

+ (2K2 -cf,-ZM, - 0 <^'WP\(Z - JV. 

In the last equation the choice of the oC 4 . 'a, depend on the 
magnitude of M 3 a~^ and the boundary data. For definiteness 
we set 

3.73 °1 “ x 

and 

3.73 a oC 2 = §M 3 (a o -M 3 )“ 1 for \>, 



or 
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3.73 b 




= M_a 
3 o 



for 



M a -1 
3 o 




Next we choose 

3 ‘ 71 * a K = 2^{l+2Mi+|M 3 (a o -M 3 )- 1 + 2K 1 2}= ^ 

or 

3.74 b K = — i~{l+2M +M a =1 +a M” X M 2 ) = K 

2 m 2 l 13o o 3 1 J 3 

whence either 

■U. 

■-} 

is given by 3*57 with K = Ky 
Subject to the remark following 3*71 it is clear 
that equation 3»70 and the proper form of 3*75 provide a 
bound for the second term on the right of 3«1&. 

By combining the results of this section, the solution 
of 3.1 is bounded at the point p 6D(y Q ) by 



3.75 a 



or 



3.75 b 



where 






6l 



3 . 75.5 



w(x„p - 7l(f)l 






+ 4)* * Xf/MT ^ 



f 




>JL 



^ (}(/( !f ■ * Y *5.)^}* + ^ { ■ IT i 7 */ 



In the last inequality the functions and ^ are given by 

2.10 and 2.39> and the constants by 3 . 42 , 3 * 62 , and 3.74* 
The computable bounds ^ ^ and ^ are given respectively 

hy 3«4®» 3*57 > and 3*75 while the function X (p) is defined 
by 3.17. 

If problem 3«1 admits a solution w which is twice 
piecewise continuously differentiable in y for (x,y)6 V, 
then theoretically the pointwise bound may be made 
arbitrarily small by improvement of the approximation 
function IP . In this regard we note that the right side 
of 3.75*5 is zero if P = L = G = 0. Furthermore, due to 
the form of the bound, it is clear that the result may be 
improved by means of the Rayleigh-Ri tz technique provided 
the problem is linear. 
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By referring to the proofs of the appropriate theorems 
of Chapter II and the methods of this section, one observes 
that the boundary value problem given in 3-1 may be 
generalized by allowing the components = a^j(x,y) to 
be non-symmetric . Such a problem was not considered since, 
even though the techniques are exactly the same, the 
individual inequality expressions are naturally much more 
involved. In the next section, where the matrix a*J is 
required to be symmetric, a method is given for symmetrizing 
a general problem. 

As a second generalization it is easy to see that the 
differentiability of the a^J, i, j - 1, 2 . . . N may be 
reduced by merely requiring that a,^ and a, ^ be continuous 
in the interior of a finite number of regions the sum of 
which is D(y) for o <y<y Q » The a^-<J are themselves con- 
tinuous in D(y) for o $y <y Q and piecewise continuous in V. 
Under this reduced differentiability, one may easily show 
that the techniques of this section are valid. 

2. The Non-Linear Dlrichlet Problem 

Our aim is to bound the solution of a normally 
parabolic problem which is such that the solution is 
initially given and thereafter satisfies a non-linear 
normally parabolic differential equation in the interior 
while assuming given values of the boundary. 

We assume that V is an N + 1 dimensional domain with 
boundary D(y 0 )USUD(o) which is such that 2. if and 2. 5 are 
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satisfied. The boundary value problem may be stated 
explicitly as follows: 



3.76 



/' 

( TAJ 



w(x,o) - g(x) 
w(x,y) = h(x,y) 



J(w) = (a ij ’ w, ± ) , . - 



: f(x,y,w,vw) (x, y) £ V 



(x,y) £ D(o) 
(x s y)£ S 



— 00 < W 7 * W 



where w is continuous in VUD(y Q ), continuously differentiable 
in V U S and is such that its second derivative in x and 
mixed derivative in x and y are continuous in the interior 
of a finite number of subregions, the sum of which is D(y) 
for every y satisfying o<y<y Q . In addition it is necessary 
to assume that 



The differentiability conditions which the solution of the 
present problem is to satisfy are obviously much stronger 
than those which were required for the solution of the 
boundary value problem considered in the first section of 
this chapter. Additional conditions must also be placed in 
the coefficients of the differential operator; in particular, 
the components of the symmetric matrix a^ J = a^J (x,y) are 
piecewise continuously differentiable in V and satisfy 2.7 
while the function 5 - 2(x,y) is piecewise continuously 
differentiable in y for (x,y)£V and satisfies 



^lim w(x,y) - g(x) 
(x,y) —>(o,y) 

(x,y) 6 V 



3.77 



lim w M (x,y) = g i (x) 

Jx,y) £ V 




i = 1, 2 . . . N. 
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3 * 7 $ min 2 (x,y) =m ? >o. 

V * 

The boundary data f, g, and g,^ are assumed to be integrable 
and square integrable over their respective domains of 
definition for any bounded continuous function w which is 
continuously differentiable in x for (x,y ) 6 V. Moreover, 
the function f satisfies a Lipschitz condition in all but 
its first two arguments. Hence there exist positive numbers 

^ /\*-f 

M, and M , i 1, 2 . . . N, such that 

3.79 |f (x,y,w 1 , v w-, _)-fXx,y,w 2 , v w 2 ) | < M, ^-w^ + M i | w x, i~ w 2, ± 

for (x,y) 6 V, - < w, vx < 00. 

We remark that the problem defined above and the additional 
problem which we shall obtain from it by introducing an 
arbitrary function are very similar to a problem for which 
Friedman, [10] and [12] , and Zeraglya, [ 27 ] and [ 28 J 
established the existence of a solution. 

The desired bound for the solution function will be 
obtained at the point p£D(y 0 ). As in the previous section 
we assume that Z satisfies 2.9 at p (notice that 2.8 is 
satisfied due to the differentiability of the a 1 ^) and that 
f is bounded in some neighborhood of p which is contained 
in V for any continuous function w which is continuously 
differentiable in x for (x,y) 6 V, that is, we suppose that 
3.5 is satisfied. To obtain the desired result we choose 
an arbitrary function - iP(x,y) which is such that J(4) 
and If (x,y) approximate f (x,y, if , 7 L ?) , g(x) and h(x,y) 
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respectively. In order that the divergence theorem may be 
applied in what follows we further require that be twice 
piecewise continuously differentiable in x and piecewise 
continuously differentiable in y for (x,y)6V. Also we 
require that the be piecewise continuously differentiable 

in y for (x,y) fc V. 

Following the techniques of the first section of this 
chapter we define 

3-80 (x, y) = w(x,y)- ^(x,y) 



and then compute 

/j( f )=f (x,y,w, v w)-f (x,y, ip , 7<f)+F(x,y) (x,y) 6 V 

3.81 < f (x,o)=g(x)- if (x,o) = G(x) (x,y)6D(o) 

\^(x,y)=h(x,y)- if (x,y) = H(x,y) (x,y) 6 S 

where 



3*82 



F(x,y) «= f (x f y, if , if ) . 



For later use we note that 

3.83 kmlSMjM'l +M l |%. | 4 |F(x,y)| . 

With J('f) and Y defined respectively by 2.6 and 2.10 

Jr 

we note that 3»l5 is valid for the present problem. Sub- 
stitution of 3 . 81 , 3.82, and 3*83 into that equation then 
yields after transposition 
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- xtf>)\ * l^w (XT ^T(Xp)^V I 

V(^) * 

f|iw ii/ {f6*, 

S >5 



X(p) = V (x 0 m 0 ) - fj" H i'fds-&w, 

c i *t*io V(n) r 

+JIy 0 ?G^.s - J T/w y„?« /s. 

to) p $ 3 p 

The volume integral on the right of 3.85 is proved to 
exist by the argument associated with equations 3«l8 and 
3„19„ Hence w(x 0 ,y 0 ) is the only unknown tern of the 
left side of 3-84 and our desired error bound is then 
clearly available as soon as the expressions on the right 
of 3° 84 are bounded. Since the unknown volume integrals 
will be bounded in terms of known quantities and the 
integral of the square of the conormal derivative of •f' 
over S, the last term on the right of 3*84 is considered 
first. For this purpose use is made of 3*8l to rewrite 
2.113 as 



3.8i| 



where 

3.85 
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3.86 



where % 
write 

3-87 



where 



+ jlf f ^{2 5^ + ft'^ 1 w j 7w)-f(i^ i «j?i(i)+f : (^-^} 5 /v 



is given by 2 . 1 l 4 « Since f satisfies 3*79 we may 



^ ^7u>) + «/v 

4. ( 5 4)^v +4 

+m«\?jv * iWF i jv 

V * V 



($5„ * I? + ?; A + T?- ttJ ^Jv 




+■ i W FVv 

V 




= max 
V 





3.88 
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and CC^ is any positive number. In the last inequality the 
constants M- and are given by 3«25 and 3.34 respectively. 
Hence 3*86 becomes 



3.89 



5 




(c+Hl 

{ L i + «.a 0 






iXT fUv. 

V 



Due to the differentiability of the solution of 3«8l, we 
may write by means of the divergence theorem 



3-90 



ffj ^ T(f)dv <-/// «.'* %Jm 

+ ^ JT/ 



Application of 3*83 to the left side of the last inequality 
yields 



3.91 



1 1JT l| TO)i v N JT/ fi, | * 1 1 fi *|^| f | f | } Jv 

< I #f(HT* + & f ^ * 1 nr($)^ 



^3. 4*" a ' ^ ^ t #jt/ ( IpV/ + h - 3 nr. 



where the i are arbitrary positive numbers. Subsequently, 
when we use 3*9^- to rewrite 3 * 90 > the arbitrary constants 
will for definiteness be defined as follows: 
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3.92 





iri2 

T 



The first term on the right of 3. 90 is bounded by noting that 



3.93 



-fff 

V °6t) * 






M B 






0(o) 



>4 J K 



H 



H 






where 




ij 



Using 2.102 the last term on the right of 3*93 may be 
rewritten and bounded by 



3.95 



- i H a \ 4/K (Hfjs 



s V 



where 

3.96 




S 




To estimate the second term on the right of 3*91 we first 
write for an arbitrary positive number (£2 
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3-97 ' ff TvJs i T JT (\vfji 



S t\i- 



and then, almost everywhere on S, define the vector 



3-98 



(T 



1 m N 

9 O * • ± $ ± ] 



where T is tangent to the surface S, is orthogonal to the 
vector (t^o), given by 2„97j> and has unit Euclidean length. 



X06# 



y~ (T 1 ) 2 + T y2 



3. 99 (T x )“ + T* = 1 

I 

We denote a directional derivative in the T direction by 

ISf - ... -ri \!t -rlr 

3.100 ^ ^ • 



Since the space is Euclidean, 



3-101 



i 

n = n. 



for every i 



whence 



3.102 



97- 






Transposing and squaring 3.102 yields 






1 
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3-103 



Since 2.5 

3.104 

and hence 

3.105 



where 

3.106 







(fi 1 ) \ 



a *% 






+ 




- trj (* 0 *l) ; (|?) i ^/K 





is satisfied on S, 

1 - n 2 > 0 

7 

from 3.103 







N } v)j. 5 4- £fr[l 

^ 5 






By applying the above results to 3«90 we have 



3.107 




In 3.107, let 

3.108 CC 2 - V~M^ 

which minimizes the coefficient of the integral of the 
square of the conormal derivative. For what follows we 
assume that Z V ) Is a strictly positive number; 

when this Is not the case, obvious modifications are made 
to obtain the same result. We set 

3.109 0( 1 = (M 9 + 2 \Tm^ 0 ) -1 



and then rewrite 3*89 a s 
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3.110 




iMy 









V 




SfJ fVv 

V 



*• iJT ^a\js 

~7. XT'” (~t)Js + JJV [^( i-fflp] 7ff ) «/5 

Mo ~k iUn, tt j , a ^^ Js ], 

3f<& 



If the last term on the right of 3.110 is neglected, then 
the only unknown terms on the right of this inequality are 
the second and third. To estimate these unknown terms we 
make the transformation of variables given by 3*20 and 
observe that 3*29 is valid for our present problem. In 
place of 3.35 we write 

3-in ( 1 - r^) JTJV*#. % Jtv 

0 y *■ t 

- T'^Z(yv)^ 4 If H Vjtp -itt 
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In the last 


inequality, let 


3*112 


a 


3*113 


S = 1 


3*114 


K = ft, * ik ^ 3 [c, + 






3*115 


= T { c , + *„*»»’*> i 3 ( Mj+lfR \^) + + 



and then substitute the resulting inequality in 3.110. 

We note that the following terms appear on the right side 
of 3*110 after this substitution 
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3-n6 - [c + (Mf 1 2 Vn,) + ^ + 






£(%) 

-A™. (a/ ? +-2V7^) 1 ff a '*?.■ f ■ Js 

W, * 

t<3. 

■•'if -^JOT 

?•»?« i(j) ^ ^ ’ C(^) 



where M^ and are defined by the above inequality and 
3.117 5L = min Z(x,y) > 0 . 

2 D(y 0 ) 

Since we are seeking an upper bound for the integral of the 
square of the conormal derivative over S, it is clear the 
unknown terms which appear in 3.116 may be neglected in the 
revised version of 3*110. In the original statement of the 
boundary value problem, the solution was assumed to be con- 
tinuously differentiable in VUS; if we make the stronger 
assumption that the continuous differentiability holds in 
VUSUD(y 0 ), then it is easy to obtain an upper bound for 
3.1l6 which is less than zero. In what follows the stronger 
assumption concerning differentiability is made since, if 
this condition does not apply, then one may clearly 
neglect 3*ll6 in order to obtain the desired result. We 



define 
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3.118 B(y 0 ) = D(y o )=D(y o ) 



and then write 



S ■f *■ n - JT f . i p*Js + a XT f Vv.ij 

B(>) ‘ OCtf U %,) " 






where 



3*120 




Next we set 



3.121 



(1 2+ CC^)a Q oC^ 





and solve for 




to obtain 



3.122 



« r = -k(-*i + 




~\iV 



)• 



Prom the above it is clear that 

3.123 ~n„ If «*V, %: Js - m., // 

mi. oftt> 1 Wtf 

*<1. 




«„*„(- + + V <t£?) 
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The inequality on the right of 3*123 is strict since in 
Chapter II the auxiliary functions f* were chosen so that 
f*n^ had a positive minimum on S. 

By combining the above, we are able to rewrite 3*110 
as follows 

3- 12 4 M$TJls*'tK*L M f'jv 

5 3 y 




■{<<* &*f> [* k , 

+ If/ F Wjo ^ 




= IS 
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where OC ^ is given by 3»ll5 and the computable constant 73 
is such that 

3.125 13-0 

when 

3.126 p = h = G = 0. 

Inequality 3*124 is a bound for the unknown factor 
of the last term on the right of 3*84* Bounds for the two 
remaining terms on the right of that inequality may now be 
obtained using the procedures of the first section of this 
chapter. With 13 known it is clear that 3*57 becomes 

3.127 JU, JJJ $ J(P)~ 

VI: ^ 

+ 1^ $ Z } Jv 

^ jJT ii ^ [v k (y 0 -i)y$ f 13 

$ 

- !f H &vf> [-- (.%- '£)] J s - Js 

The right side of 3*127 is computable and hence the desired 
bounds may be obtained in exactly the same manner as was 
employed in section 1 of this chapter. The results of these 



78 



where <jC ^ is given by 3.115 and the computable constant IJ 
is such that 

3.125 18 - 0 

when 

3.126 P •- H = G = 0. 

Inequality 3*124 is a bound for the unknown factor 
of the last term on the right of 3 . 84 . Bounds for the two 
remaining terms on the right of that inequality may now be 
obtained using the procedures of the first section of this 
chapter. With 13 known it is clear that 3*57 becomes 

3.137 itm Iff {if J(*P) 

+ Z % % J/ 

? JJ'V 7 H t13 

s 

- [f V5 H - (( i H [?_ kfy-fljjs 

1l3 rz ' &2 ^( ih ^- cls = } 

The right side of 3*127 is computable and hence the desired 
bounds may be obtained in exactly the same manner as was 
employed in section 1 of this chapter. The results of these 
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operations may be written as follows: 

3.128 

< { xr;[-jc?’)]' , j(^v}Vf sss[mT(?FF) x . 

<- {$*. 

+ iix t [2 \ 

4{h^Js}Ht} k 

where the functions V and are given by 2.10 and 2.39? 

P 

the constants K^, and K by 3.114, 3*62 and 3.74, 

. J /ys 

equation 3*127 defines ^ ^ and the constant $ is given by 
3.75. The computable bound "fS is obtained from 3*124 while 
X (p) is defined by 3*85* As in the problem of the previous 
section, we remark that if the function f is such that 

3*129 f(x,y,w, vw) = f (x,y) 

then the result lends itself to improvement by means of the 
Rayleigh-Ritz technique. 

In obtaining the conclusions of this section, the 
role of inequality 2.113 was basic. By referring to the 
derivation of that inequality it is clear that the symmetry 
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of the components of the matrix a*3 was essential; hence the 
consideration of a method for the syrametrization of an 
arbitrary problem Is appropriate. The result Is obtained by 
defining 

3.130 a iJ = ■|(a i 3 + a^ 1 ) 



and then, since 



3.131 aiJ = a 1 ^ f, J± 

almost everywhere for (j,y)6V, we may write 



3.132 



= <1 K * « 4 % 



where 3*132 Is also valid almost everywhere in V. Due to 
the foregoing, the differential equation of the boundary 
value problem may be replaced almost everywhere by 




= f ')• 



3.133 
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Prom 3.133 it is clear that £ satisfies the same conditions 
as £ and thus the original problem has been symmetrized as 
required. 

Before concluding our consideration of normally 
parabolic problems, we remark that when the differential 
operator has the form 

3.134 J(^) = (a ij , 

the transformation given by 3*20 yields 

3.135 J(?) = (a lj 9^) ± - (^ + K )•? • 

Hence, if the minimum of the function ^ over V is a strictly 
positive number, then the constants of either this section 
or section one are reduced by the amount of this minimum. 

In order to appreciate the improvement which results under 
such conditions, the reader is referred to the example 
presented in Chapter VI. 



CHAPTER IV 



DEGENERATE PROBLEMS 



1. The Dirlchle t Problem at a Parabolic Point 

The problems which have been considered thus far have 
been normally parabolic In character. We now turn our 
attention to problems which are normally parabolic at certain 
points and degenerate at others. As the above heading 
implies, this section pertains to the derivation of a bound 
for the solution of a degenerate boundary value problem at 
a normally parabolic point. We seek to bound the solution 
of 



where w is continuously differentiable for (x, y) € V U S U D(y Q ) 
and is such that its second derivative in x and mixed 
derivative in x and y are continuous in the interior of a 
finite number of subregions, the sum of which is D(y) for 
every y satisfying o<y<y C) . In addition it is necessary 
to assume that 



<W 

\W 




(x,y) 6 V 
(x,y) e D(o) 
(x,y) 6 S 



4.1 



W ( x , o ) = g (x) 
W(x,y) = h(x,y) 



4*2 



lim w'x,y) = g(x) 
(x,y)— 5>(x,o) 

(x,y) e V 



lim w, ± (x,y) - g^^ (x) 
(x, y) — ^(x, o) 

(x,y) 6 V 



i * 1, 2 ... N 
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The components of the symmetric matrix a*J satisfy the 
conditions imposed in the second section of the previous 
chapter. The functions q and 2 are piecewise continuously 
differentiable in y for (x,y) 6 V and Z satisfies 

4*3 Z(x,y)>o (x,y)6V. 



Prior to stating the additional conditions which are to be 
satisfied by the functions Z and q, it is necessary to 
define a subregion E of V by 



44 



E 



= {(x,y) Z(x,y) -0, (x,y)6v}, 



We suppose that q has the form 



k-5 9 = Z I t I 

where the constants £ 




i = 1, 2 . 




satisfy 



d r 

4.6 I £ r > , 

I~-l 

and the functions q^., i = 1, 2 . . . , fulfill certain 

conditions which will be prescribed in what follows. Let 
the coordinates x^, 1 - 1, 2 . . . N, be chosen to satisfy 

4*7 min max / (x^-x^)^ = R 2 

(x»y) (x,y)6S 1 

where o i.yi y Q . For N = 1 we suppose that 

4.8 a m = min(o;q> -a o (4R 2 ) 

E S 
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with E C 7 defined by 

O 



i|.9 E^=|(x 5 ,y) ||(Xj,y) -p || < $ for some p <= E with S>Oj/ r |V 



In 4*9 the symbol || || denotes the Euclidean distance in 
the 2 (and later N + 1) dimensional space of the problem* 
If on the other hand N = 2, then q must satisfy 



4.8 b 
where 
4*10 
and 



m =minio: 



E, 






r i*i> 



C4h x ) 1 ] 



1 = 1, 



A 




1 

2 



4.11 R— “ max r . 

1 V 1 



Finally for N •? 3 



4.8 c 




l 



Since the techniques which are used to obtain solution 
bounds for N = 1 or N = 2 are very similar to those which 
are employed for N$'3» the method shall be presented in 
complete detail only for the latter case. Thus for N$3 
we define m by 



4.12 



0 < -4 E I m I 




m < 1 



If N = 1 or N = 2, similar constants are introduced. 

Previously we required that Z be piecewise continu- 
ously differentiable in y for (xj,y) 6 V; in addition to this, 
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(x,y)€E where E is the boundary of E, and f^-o on Eflfl-V). 

In 4*1 the boundary function f » f (x,y) is continuous 
in y for (x,y)6v and the first derivative of f with respect 
to y is continuous in the interior of a finite number of 
subintervals, the sum of which is the line (x, y) for every 

— — "N f 

fixed x where (x,y)6V. In addition we suppose that ~-f- 



Is integrable and square lntegrable over V and that g and 
g,^ are Integrable and square lntegrable over D(o). By 
definition it Is clear that f Is lntegrable and square 
Integrable over D(o), D(y Q ), S and V. 

It is now our aim to determine a bound for the 
solution of 4.1. As the heading of this sootion indicates, 
we assume that 



and that 2.9 is satisfied at the point p 6 D ( y 0 ) at which a 
bound is desired. The arbitrary function ^ (x,y) Is 

chosen as In the previous chapter and the function f (x,y) 
Is Introduced where 



4-13 



Z( x,y) > o 



(x,y) * D(y 0 ) 



4.14 



l /'(x,y) ~ w (x,y) - if (x,y) 



We compute 



4.15 




JCH * P(x,y) 

Y (x,o) = G(x) 
if (x»y) ~ H(x,y) 



(x,y) 6 V 
(x,y)£ D(o) 
(x,y ) £ S 



Let 
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4-l6 J( t / / ) = (a 1 ^), j - q r + -^L (Z f), 

and then with 6 given by 2.10 we note that 3.15 is valid 

P 

for the present problem. Prom this it follows that 
^• 17 

ffs w T(&>) Jtv\ 

n* v(^ 1 

where 

^° 18 %(p) ~ LPfrojfy) ~ SfH vf ds 



l.'j OpFjv + JT YpiGdS ~ fj/n . 
a -*8. r 0 (o) r 5 > p 



Since P is bounded over V, it is clear that the volume 
integral in Ii.,1.8 exists? thus the desired bound is available 
when the terms on the right of 4 .17 are estimated. We shall 
first consider the integral of the square of the conormal 
derivative over S and then use this bound to estimate the 
unknown volume integral. For this purpose we make the 
transformation of variables given by 3*20. The boundary 
value problem then becomes 



8? 



4-19 



?(£)= 5 *?[*(*-*)] (-,^ev 

?(.V) » S -t^o (Kg,) (^)C-D(oj 

y (^3) = « [k (^-g)] (^S 



where 



4.20 



q — K 2 + q , 



Dae to the fact that .Iz 

is continuous in V for (x,y)€ E 
0 y 

and 4.8 is satisfied, we may choose K > o so that 



4-21 c 
for N ^ 3 and 



/ KH . - /WWW 

I y 



IHV 






1 1 > a. _ a N 

4.22 min (4 kZ - x 4~) > - (1-m) -2-_ 

V 2 5f 

where 



4-23 



q = KZ+ a I q T = 2_ £ x (-tos+q ) = £ 

1 1--1 & 6 1 1 ] 






In the above we have made use of the fact that -2-= = 0 

bf 

on E. Similarly the constant K>o is chosen so that 
4.21 a m = min{q-^2K > -a o (4R 2 )“ 1 } o] for N = 1 



m = min|q-^2K > 



and 






I 
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4 . 21 b Aj- - { o ' \ ( % l - ') > -^(V^ r ) -/ ] 

I 1 Ijij., ,(J. Oy/^d 7, 

With K>o as chosen above we have from 2.113 



b-zk < ss(%YJ: 

C * / 






+ ^' v JTJ ^ XCT 

+ ^3 ^ ^ ^ V' 



where 


8i ^ is an arbitrary positive number and 






r> tL ~) 


4-25 


m 13 = “ x 1 


1 z iff 1 ) 2 } 


4.26 


M, 1 - max ■* 

x 4 v 


l5 2 } 



In order to estimate the unknown terms on the right of 4*^4 
we observe that for any set of N functions g^, i = 1, 2 . . 
. N, which are piecewise continuously differentiable in x 
for (x,y ) 6: V, it follows that 

4* 2 7 " -S !?,*// -* 2 jXT^ ^ f jv. 

By setting 



4.28 



g i = (x^x 1 ) 
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In 4*27 » it is easily shown that 



4-29 




Substitution of this inequality Into 4-21+ yields 

4-30 UiUfds 




t-S 3 JJf ^ JJ/? (ff ) dv. 



To bound the unknown terras on the right of 4*30 we write 

4-31 f(V)Jv = JJ 5 F Ji - Xf OFo^p [ZK^Js 

* 0(1,) Df) 1 0 






. ^ 

+ ^ 3 HF^o^KC^-^yS «- J3T 
+lf l| 5v^ -TEjT £ ^ 4 
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+t 



XT 




Due to the assumptions which were made concerning the 
differentiability of the solution of 4° 1» it is clear that 
3.93* 3-95* 3-97* and 3»105 are valid for the problem 
presently under investigation. Substitution of these 
results into a transposed form of 4*31 gives after making 
use of the fact that ‘A = V on D(y_) 

■ip* 



4‘32 



f +JJ ( 

v d * %) U 



^zK^l/7 T /d ) 

$ 



+ Sf F l Js 
J>W 




. vr 







4 X ^ 

* 



dv + % 



where we assume that (M^ + 2 Vm^q) is a strictly positive 
number. If this is not the case, obvious modifications in 
what follows will yield the same result. In 4-32 it is 



clear that 
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4*33 



+ t///[r^ “Kf]^ [2kc^-^)] ^ 



f 4- ^ ^ & P ^f( ZK 1io)^ ~k H l Zyf[2><(y u -y)]M 
** * ^ ^ ^ + +A 0^ (* V)/* , H 



“iXf^ i* ( 3t) >v^0 







and 

U*3U M 1S = ““ |"-|| j 

We now ahow that due to the form of q and henoe q it is 
possible to choose so that tho terms In the braces on 
the left of 4*32 may bo bounded from below by a computable 
oonstant. On the surface D(y 0 ) we have by means of the 
divergence theorem 
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1».3S 






+ 2 // 



where Bfy^) is the boundary of D(y Q ) . For N 2 3 the form of 

4.35 leads us to define 

4.36 g 1 = 

and then for every I obtain from 4 '35 

4-37 || ^ -* j 5 

5 ( 3 .) 1 



^ If ^ ^ 'H h 



where cC is an arbitrary positive number satisfying 



4-38 



ot > 



1 

N-2 



We choose cC so that the coefficient of the last term on 
the right of 4*37 is minimized. The above equation thus 



becomes 
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4 - 38 c 



// 

®fw 



4 (Zjfi, SS a. '*({ V' ofa . 
‘ 8(4) 1 “ (W2 K Pfo) * V 



For later use we observe that It is also possible to show that 



4*39 o ffj P^Jtv 



£ l-i $£f*\ *- 



For N = 1 we 


have 








4.38 a 


S Ms ( 


X (*-*) </ 


+ 


vi 2 r " 

-f j? 








X* 


6 w 


and 










4.39 a 


fX?Vv - 

V 


2 / (* 1^#/ 5 +• 


^ UVV t* 


Finally for N 


- 2 we may 


write 






4.39 b 


JT vV'Js 
w» 1 


w %) 


and 










4.39 b 


/// 







o v 









i 




\ 



I 
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We are now in a position to obtain a lower bound for the 
terms in the braces on the left of 4*32. The details shall 
be worked out only for N^3» Similar results easily follow 
in the other two cases. By means of the above and 4*21 c 

4-4° ff r qtjs z JTX e x ^ r *\ 2 Js Hi Ji 

m Iz/ 1 o&) 



* 



l 

I--\ 



f r - 

O /vyi 



( k/ - z ) )/if ~/vn ff ^ ^ . y., Js 

S 1 " DQJ ' * 



We notice that the line integral on the right of 4*4® 
actually a volume integral over an (N-l) dimensional 
domain; and, since N £ 3, it is clear that the integral 
exists. In 4*32 we choose 

- N 2 a 0 

4.4l £ 1 = >° 

where the strict inequality on the right follows from 4*12. 
Substitution of 4‘4l into 4*32 produces the inequality 

l *- lt2 Iff i(s%) fM„)e*pfekty)lf (MfJi 



+ZK‘ 



A 



(/--) to 



+(&***' 



'» N% 
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Cv-2^ y v 
%) * 












f (*'-/) /n. H ! <Ji f A 

S(g») 



In deriving the above, an inequality analogous to 4*^9 over 
D(y 0 ) is used to bound the integral of the square of ^ over 

D(y 0 >- 

So that 4*30 may yield a non-trivial bound for the 
integral of the square of the conormal derivative, we choose 
0C 1 as follows? 

4-43 = |m 3 (M 9 +2Vir o )“ 1 exp (-2Ky Q ) 

Inequality 4*30 then gives 



4-44 



jT(s>) Js i ^) 3 i c , +■ *■ + **{>(**$,) 

= % jyjv'e. t 1 



where 




: 






I' 

I 

i 









96 



4-45 



^ ^ + ^yj H H H lj ^p [ 2 K O30 ’ ^] ^ ' 



+ ^ -tt" f t l^tZVM,,) bp (-Z^{^ l [u\ (l-wj\ ' 

•iT^Vi + ^ - >-£ r ™ r 6^2 
+0-^) ^ r f/^ij 

— o(ij) 

The unknown term appearing on the right of 4*44 * s estimated 
by using the divergence theorem to write 



Stf i*z A it = - ssi ? T(p)dv 

v r v 

+ (f ? Iv ^ - JXff $ 7 iv ' xXf * 

-I ifr Li ? Vs * 4 JX i Ms t+SffM Mv 

^ i D(o) V 4 

% '(‘'t -f/ ITOc^x^ // T ^p[z/r f 3£»-^)]W ^ 

+2l, Iff pVp[zk^)J^ + | 2 ‘I/( 



it H z **f[ L iK(q 0 -'tt)\js 
s 

£ 2. \ c 

~Y, ^ 

+"* iTJ 0. * h + x (z K ^ 
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- 4 l/-^ 2 //Wy^K + i: (j-^) £fJV ^^jjv 

where 6 and £ are arbitrary positive numbers which shall 
be suitably chosen in v/hat follows. To obtain the above 
it is clear that 4*21 c, 1^.22, 4-29 and 4«38 c have been 
used as necessary. Inequality 4*4& now transposed so 
that the unknown term on the right of 4«44 m& 7 be bounded 

Iff a'* % £ Jv 

v K i 

*kSf( I v) Z Js f TT^HWp [n(w)]Js 

5 5 

+ ( t r L ) f° V-2*) *&,- *)] J S 

2 JXTf jz/c(^ 0 -^)J Jl/ 

’^r feI )\\- H Js 

+ 4 - T 1/^* 

+ (/-w) TP // (x l - x‘) XI, [2 K (^ a -a)j 5 

iS 

= T-^(lv)^S t^ y 









'I 

I 
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where ^if Is defined by the above equation and it is clear 
that we have set 

448 6 1 « (1^2 !) 

Since m<l due to 4.12, inequality 447 provides a non- 
trivial bound for the Integral of ^ By comparing 

444 and 447 we are led to choose 

449 £ 2 » 

then 4 44 becomes 

4-50 [J(Wjs i t\ f %*(i 

+ = “6 

v/her © f3 la computable and la such that 

441 13-0 

if 

442 P - H * 0 - 0 

The above la a bound for the unknown factor of the 
second term on the right of 4* 17* Wo next bound th© volume 
Integral on the right of that Inequality in terms of 
computable constants and 440. Let 

4-53 A 1 /') » (aU ^ (1 ) - 2q^ + ^ (BV) 



I 






{ 
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and consider the function "f 7 given by 2.74* Since the lower 
order terms of and "j are identical, we may deduce from 
Theorem 2.3 that the arbitrary constants of 2.74 ma y b© chosen 
so that 

4*54 " ^ M and £0 for (x,y)£V 

and 

4.55 o 

for every (x,y)€V such, that o<£< Y 0 ~Y£J 0 ‘ Prom the 
proof of Theorem 2.1 and the definition of f and 2., it 
follows that 

4.56 

By means of the divergence theorem and the differentiability 
of we have 

b-57 J2m. ^(T)} 4V-SIU 

- SI fa ^ Tv } ds - ^ z Js 

S. s J 



■f 



r/^2 fVi 

Wo) 
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and by expansion of the last equation the following inequality 
is established 

4*58 iUv IJJ h « -zjjft'afy Jv 

V(t) ' A V 

h^o 

J-fffr 2 - Vi H z J.s t G l Js + 13 

S d D(o) 

< ZM 0 1 SSifF z - ^t-/y) Z^} H z Js +ffV-^Vi-/-$- 
s y m o 

where £) is obtained from 4*4? after application of 4»50 
and M is given by 4*54- By an argument similar to that 
which was used to establish 3.64* we are able to write 

4-59 ^ 

and finally from equations analogous to 3*86, 3 » 87 » 3»88 
and 3*89 we have 

4* 80 iff 

+ {£ f Slf[-T*^)J^^ z Jv} k ). 



The existence of the integrals on the right of 4«8o is 
established by a discussion similar to that of the previous 
chapter. 



Inequalities 4 ->51 and 4*^0 give computable bounds 
for the unknown terms on the right of 4*17 and hence the 
desired point-wise bound is obtained from the inequality 
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4.6l 



I w k Ja „) - { JTJT- J W' ■ 



+fj 



where the functions X and t 7 are given by 2.10 and 2.74 

.H 

and the computable constants ~i3 and ^ are obtained from 
4.50 and 4 «5>8. As in the previous problems, we observe 
that the right side of 4«&1 is identically equal to zero if 
F = G = H = 0. 

During the investigation of the problem of this 

section, an attempt was made to generalize the method by 

^ -3 

omitting the requirement that be continuous in V for 
(x,y) 6 E. Such a generalization has not as yet been achieved 
however, if the domain V is restricted as indicated below, 
then the desired bound is available. 

We suppose that the given boundary value problem is 
identical to that stated at the beginning of this except 

o » 

that need not be continuous in V for (x,y)eE. In 

ay 

addition, however, we suppose that V is a cylinder (that 
is, D(o) = D(y) for o^yly^). It is evident that equations 



I 



i 



( 
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4-l4 4*20 are valid. We choose K>o so that 



4.62 



m =minfo s r^q > ~ra . (N-2) 



V 






I^I 






I = 1, 2 



J) 

• • • 



for 3 where it is clear that K>o may also be similarly 
defined for N = 1 and N = 2. In what follows we shall 
consider in detail only the case N£ 3 since similar 
techniques would give the same results in the other two 
cases . 

The integral of the square of the conormal derivative 
of ^ is obtained by botinding the unknown expressions on 
the right of 4*30 in terms of the given boundary data. Prom 
the divergence theorem and the differentiability of and 
the coefficients of the differential operator, we have for 
an arbitrary constant "a” 

4,63 JjfT% bUri)]^ 



' ■ if 6F **P [( wk )^o]ds 



Hit*) Iff 4*f> [(«**) (ft-d] A' 
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*(“ * f "2 If 

±h, Ifft^'H ^p[^(a+2k)C^-^)]^ - -kXf ^ ? ds 
5 ^0 

+ H4j ' *$Tf ^[*w]^ 

^JJTff fifW^ - 



The above inequality is transposed and use is made of l±.2 9.< 
4*39 c> 4-^2 and Schwarz's inequality to obtain 

su\ *> + f 

S T^( 5v)^ + 4 lf[t|~ Krt ] 



^ frW ' 






•f 



i;. ST F ((i 4 J/ fpi%^p[(a(-2K)Mj^ 

«*) »0 1 d J 

( A4 ' K ) F ] 



In 4*30 and the last inequality, we set 



4-65 


i! 

M Lf 


4.66 


<^1 - 1 




N 2 a Q 


4-67 


£? ~ 2 ■ (1 " m) 

2 4r 


and 




4.68 


/ __vir M 8 

a. -(I-#*) ia t + hiq o + 






Substitution of the above into 4*30 yields 

k ' b9 in I tyjs -2% 



’S 






^k o 0-^) // ^ 



where a is given by 4*68 and K>o is chosen so that 4»^2 
is satisfied. Since t3 is a computable constant, the 
desired bound for the integral of the square of the conormal 
derivative has been obtained. A bound for the unknown 
volume integral in 4*17 is again given by 4 * 60 , and thus the 



io5 



desired pointwise bound is available. 

2. The Dirl chle t Problem at a Degenerate Point 

We seek to obtain a bound for the solution of 4°1 a t 
a degenerate point; that is, at a point p which is such that 



At a degenerate point the author was unable to determine a 



normally parabolic point. Because of this reason, additional 
conditions had to be placed on the boundary data, the shape 
of the domain and the region of validity of the differential 
operator. 

Let the problem for which a bound for the solution is 
desired be that which is proposed at the beginning of the 
last section (specifically, from the heading to the paragraph 
immediately preceding equation 4 « 13 ) with 4*70 replacing 4 * 13 * 
In addition we suppose that V is a cylindrical domain, 
is twice piecewise continuously differentiable in y, and the 
differential equation is satisfied on D(o) . Moreover we 
assume that the derivative with respect to y of the conormal 
derivative of w is continuous in the interior of a finite 
number of regions the sum of which is S and that the fol~ 
lowing are satisfied 



4-70 



2(p) = 0 



parametrix with properties analogous to those of ^ at a 

A 1 ' 





p,eo(c) p z e s 
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4-72 



Z(x,o) > o 



4-73 




£ -»(x,o) 
p 6 V 



and 



4-74 



Z(x,y 0 ) = o 



for some open set D'‘CD(y o ) for which p is an interior point. 

To shorten the argument, we consider only the case 
N £3 since for N - 1 or N - 2 similar results are easily 
obtained. As in the previous sections, we obtain the 
desired bound by introducing a function = Ifi (x, y) with 
the same differentiability as w and then with ^ given by 
4.14 we compute the system 4»l5* Next we define the function 



V given by 3*20 and compute 4*19* The constant K ^o is 
chosen so that 4*21 c and 



are satisfied. The desired bound is now obtained by using 



4*75 




( 6 . 8 ) 




to write 
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where 

4-77 d(f) = (YVJ 






and Pp is given by (4*2) of [2lJ . The unknown terms on 
the right of 4*75 are the first, second and fourth. For 
the fourth term consider that 

4-78 f r f j m}? {§ (.VUfJs }' 1 . 

Kl,) %) %) 

Using Hormander's Identity and the same techniques as were 
used to obtain 2.113, it is clear that 



4-79 


§ {.Wu s % 

Bfo) 






where 




4 • 80 


4 § 




+ 1 f n 1 f ( 

3 Btti * 


4.81 


-1 „k 

m = mm n f n. 
3 B(y o ) k 


and f 1 


= f^(x) are any set of 






V 



A 
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differentiable functions which are such that f^n^ is 
bounded and has a positive minimum on B(y^). Prom 4»32 
after application of l\.$0 and 4*47 or from 4»&4 with the 
unknown term on the right bounded by 4*^9» it is evident 
that one may obtain 

4.82 ft < k >r %. *» «/s s D 

* * 

where is a computable constant. Inequality 1^.82 bounds 
the last term on the right of 4*79> for the second term we 
write 



4-83 



D(tf * * 

~ z JT « vds + £ // ^ ds 

^ J>($ * * 

+ zSf 

Pl») * 





f XT 2 (l|) <Js - 


where 




4.84 


M-, - max ^ ^-(f 1 ) 

D(y 0 > ' 


4-85 


TL 0 - max {s )_(f 
D(7 0 ) > 



i>fo) 



SS F*/* 
0(8.) 



and 
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4.86 




By using an inequality over D(y Q ) which is analogous to 4*^9^ 
we may rewrite 4 4 79 as 



4.87 



%>) 








if: , gj 

% 4> 




3 



+4 // F//3 + 4 u 

3 ^,) 3 m * 

• ' H *4 &• 



For the last term on the right we use the divergence theorem 
as follows: 



4.88 




Since 
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4.89 







/ 







It follows 



4.90 



that 

\o.'~ 




Prom 4*4-7 after application of 4*^0 or from 4*^4 with the 
unknown term bounded by 4*^9» one has 



4*91 jJjV 1 ^ tjjv 4 

By transposing 4*88, and integrating by parts, we obtain 




+ 7v l -|X (* F*# + 



where 



iXT [-KF+ l|j </v t ^ 



Ill 



4-93 

1 1 *'" ~ 3 ^ 4 £ *&•#)] ^ f i tS 

- XJ £ J ^ r (tpt) //(* - >£ ) ^ r V {-/t// + 1^}^ f ^r^pfs 

s 

+ H/ 2 ''f ( a ‘\),j -t fr - ^} 2 ‘tffafr) ds 
+(/- ^ j ^ri Sf(fi A - i") 4 !, i~ K ^* ^ 1 '£)] ds. 

In 4*92 we choose the strictly positive numbers oC^ and oc ^ 
so that the coefficient of the second term on the left is 
zero. By comparing 4*87 and 4*92 we are led to define 

4-94 ^3 = t S 3 

By fixing the arbitrary positive constants as indicated 
above, it is clear that one easily obtains 

4.95 f (uyj£ < ^ 

_ 

where |3 is a computable constant. Inequalities 4*95 and 
4.78 provide the desired estimate for the last term on the 
right of 4*75* The second term on the right of that 
equation is bounded by using exactly the same procedures as 
Payne and Weinberger employed in [2lJ . Let^(p) be defined 
by (4* 18 ) of [2l] t it is then clear that /O and Jb- 



may be chosen so that 



4.96 

and 

4-97 

where 

4.98 

Prom the 

4-99 



and thu3 

4.100 



i M 



for (x s y)6D(y 0 ) 






for (x,y)fD(y o ) 



divergence theorem and the definition of 

If -v>d*&)Us - ff{ *i?*\ 

= - ff r z a*fr)ds + zSfv' , faM 

zj'TYJt t i3 -f »* 



{Jf-v'V'c^Wsj^ 

%) J 

+{ UH*&)r*F' a(vu s 
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To bound the unknown terra in the above expression we make 
use of 4*82 as follows 



4.101 


(Sfjs 5 ei -- % 

W 0Cio) a ° 7 


and then 


obtain 


4*102 


%) 




jo) ^ 0 ) 






where 




4*103 


M , = max { JLa* )J " V 2 S j 

15 D(y 0 )-D" 


4.104 


NI = max LgK 2 } 

17 D (y 0 )-i* J 


and from 4*92 after application of 4*95 


4.105 


if? <- 1 

* 


By means 


of Schwarz's inequality 


4.106 


is; ^ atr p )Js\ 

i {IT- * VflO}* { XT [- a ( r,fJs f l 

%) 0 (*) f 



In the above the first factor on the right is bounded by 
4*100 and 4*102 while the second may be computed due to 



(4*1 b) of [ 21 ] . To estimate the first term on the right of 
4.75 we consider 



4*107 


fj Pm) is 

~ 1 / ^ ^ ^ F ^ V 


where 




4.108 


M. 0 = max {g r 2 1 
18 D(y 0 )-D^ P J ' 



Transposing 4*75 and taking the absolute value yields 



4.109 


Jf Fff-Js -§ H | 

Ot*J r B(fr) 

+ZM ls % s +\Sf VQ(P f )J s \ 

+ lf Pf SvJM 
Shi) r 



where 4*106 and 4*78 bound the last two expressions in terms 
of computable constants. 

It is the belief of the author that the conditions 
imposed on the boundary value problem of this section are 
excessive. In future investigations, efforts shall be 
made to produce pointwise bounds at a degenerate point 
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without imposing the rather restrictive conditions which 
were necessary in order that the above techniques yield the 
desired result. 



3. The Mixed Problem at a Parabolic Point 

Thus far in this chapter we have always required that 
the matrix be positive definite and that g be a non- 
negative function. Now we suppose that 2 is strictly positive 
but reduce the conditions on the matrix by requiring that 
be positive semi-definite. 

Let V be an (N + 1) dimensional domain and let the 
notation of the problem be as indicated in equations 2.1 to 
2.4- In place of condition 2.5 we now require that 

4.110 min riy = m^> o 

We shall derive pointwise bounds for the solution of 



4.111 



< 



^ (*) 



V 

where clearly it is necessary that 



4.112 l (x,y) = 0 

for (xjyJ^S^ with 



(x,y) 6 V 
(x,y) € D(o) 

(x,y)€ S 



116 



4-113 



S = { (x,y) (x,y) e S 



f- n 2 

Vi 



’} 



The solution function w, the coefficients of the differential 
operator, and the boundary data have the same differenti- 
ability and integrability as the analogous functions of the 
mixed problem of Chapter III. As in the previous problem, the 
function f satisfies a Lipschitz condition In w, that is, 

4.114 |f U^y^J-f (x,y,w 2 ) 1^ |w 1 - w 2 | , 



As the title implies, the solution is to be bounded at a 
parabolic point p6D(y 0 ). Thus a x J Is positive definite so 
that a^j(p) exists. In addition to the above it is 
necessary to assume that conditions 2.8, 2.9, and 3*5 are 
satisfied at p. 

We use the notation and techniques of the previous 
mixed problem to obtain the following equation which Is 
analogous to 3-1& 

4.115 SfS ^ T (^ p ) Jlv \ 

W* %) * 

0 w 

t { ir Ms} 1 * \ * *$~jls } 

s s 0 



where 
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4.116 Z(p) * W f JJ L V Ji 

+ 11 Y,, a 6 is - /// r« F Jv. 

D(o) W'tf) 

The right side of 4»ll6 is computable by the arguments 
associated with equations 3«l8 and 3 « 19 > and hence the 
desired bound is available when the unknown expressions on 
the right of have been estimated. To estimate the 

boundary integral the change of variables given by 3*20 to 
3.22 is made and then by means of the divergence theorem it 
is clear that 

^* 117 J $ - 2 (/** /»? ) ; { X SI $ Js 

+ i J L z Mf p. +• (-/c^ +M l +t[tjn l/ )Jffp 2 J\/ 

S v' 

^ XT/f 1 ^ [l K 1/ + i SJ* (r Z Mf (2 «fy) M] 

where 

4.118 min 2 = nu > o 

S * * 

In order to obtain a non-trivial bound for the left side 
of 4*11? we se ^ 

4.119 \ 



and 
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4.i2o K = +-i^H v ) = A / 

Dae to the above we may write 

4.!2 1 || i V L Z ^D [ 2 ^i(^ 0 -^)] Js 

fLC/m^y' fjf 

f 2. C ; ^ K > ^ 

« tf 



In order to bound the volume integrals of 4»H5> we 
introduce the function given by 2.39« Even though the 
differential operator of this section does not satisfy the 
conditions of Theorem 2.2, an inspection of the proof 
reveals that the conclusions of the theorem are valid for 
the present problem. First we observe that constants 
0 , , and associated with equations 2.36 and 

2.37 exist since J is normally parabolic at the point p 
and thus a*J is positive definite in a small neighborhood 
of p due to 2.8. Because 2 is strictly positive in V 
it easily follows by remarks similar to those which were 
employed in the proof of Theorem 2.2 that the constants 



, 













oC and Sr of 2.39 may be chosen so that 



4.122 



J(^ ) < o, <fio 



cwT . 



(x,y) 6 V 



and 



4.123 J(^ ) < <0 , (x,y]6{x,y (x,y)6V,o < £ < y Q 



Due to the foregoing, equations 3«49 to 3*75 are 
valid for the present problem with 



4.124 



M 2 = m 1 = X) . 



The bounds for the unknown volume integrals of 4*H5 are 
hence available by means of the appropriate equations of 
Chapter III. After application of these results one 
obtains 

1j. 12£ - Z(jo)| 

< { SSH-KT)]* ffyflv}'- 

■(Mft-r &)]"(? Ff+p h^-$}h} h 

^ ^ + ty ( QpF ) x ^f i)]^T) 

tf XT ( Sv* + ^ { 13} ^ 
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where ^ p is given by 2.10 and Is defined by 2.39 subject 
to the remarks of this section. Using the results of 
Chapter III we obtain that Kp and K-^ are given by 3*62 and 
3.74a while the computable constants ^ and are defined 
by 3.57 and 3*75a. The bound $ is determined by the methods 
of this section and is given explicitly by 4«121. 



CHAPTER V 



THE ELLIPTIC PROBLEM 



1 . The Neumann Problem 

In [ 21 } , Payne and Weinberger obtained pointwise 
bounds for the solution of second order interior elliptic 
differential equations with Dirichlet or mixed boundary 
data. The occurrence of the zero eigenvalue of the free 
membrane problem prevented the authors from extending their 
results to Neumann boundary value problems. By following 
closely the procedures which were used to obtain a bound 
for the solution of a normally parabolic mixed problems,,, 
we obtain such bounds for an elliptic problem of the same 
general type. 

In what follows we shall use the notation and many 
of the results of section 6 of [2lJ . We seek pointwise 
bounds for the solution of 



5-1 



B (w) = (a ij ”w .) 



-b'Sr 






■qw = f (x) 



x 6 L 







x 6 B 



where D is a domain of Euclidean N - space with boundary B. 

In order that the divergence theorem be applicable it is 
necessary to assume that w is continuous in D, continuously 
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differentiable in D and is such that its second derivative 
is continuous in the interior of a finite number of 
subregion, the sum of which is D. In addition it is 
necessary to assume that the solution of the boundary value 
problem is such that 

5.2 JLim a 1 J '(p)w (p) = a^ (x) w , . (x) i= 1, 2 . . . N 

£ ->x £ B J J 

P £ D 

o t A Jt 

We suppose that the components of the matrix a 1 *' and the Jjr 1 
are piecewise continuously differentiable in D while q(x) 
is bounded and strictly positive in the same domain. As in 
( 6 . 3 ) of [ 2 l] it is necessary that 

for all real numbers ( § Q • • • where equality holds on 

the left if and only if all the numbers are zero. Let 

5.4 min q(x) - m, > o 

D 

where m^ is positive since q(x) is strictly positive by 
assumption. Prom ( 6 . 8 ) of [2l] we have 

w(p) - f(p) - JT { ^ B (w-4) -(w-ip)B (/Jo)} A 

n )~ PA(^)}4S 

8 r 



5.5 
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where 

5.6 B< Pp) = (aij r .) . + c a 1 r ) -q Pp 

P» J » i p , i 

with Pp given by (4*2) of [2lJ and $ defined as in the 
second paragraph of section 4 of that reference. We apply 
Schwarz's inequality to the first boundary term in 5*5 an< 3 
obtain 




where, as in the previous parts of this thesis, we have set 



5.8 f 

Let f^(x) . . . f N (x) be the set of functions introduced in 
sections 2 and 3 of the work by Payne and Weinberger, that 

-f 

is, the set f (x) is piecewise continuously differentiable 
in D and is such that f^n^. is bounded and has a positive 
minimum on B. For later use, we set 

5.9 nip = min f^n. 

B K 

and then bound the unknown factor on the right of 5*7 by 
writing 

= k x XT %- Vf'Jv t 4 Xf Jv 

1 D ^ ' V 1 



5.io 
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where 




and £ 1 ia an arbitrary positive number. By means of the 
boundary value problem and the divergence theorem. It is 
evident that 



S-13 



JffeJv = - JTv'/V, A + S Js 
£> ' 0 4 * <? 

D £> * 



where 

5.14 F(x) = f (x)-B( </>) 
and 

5.15 l(x) = -?(x) - . 

By means of 5>*3 and a transposed form of £.13, the 
following holds for an arbitrary positive numbers £ ^ 
and £ ^ 
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5.16 



A XT £51 f 

' If &* %.%■ h -h If& A %■ pdf' +Sf(e- %) V z dv 

0 * * D *■ D 0 A/ 




IS FtoVv 

D 




+ 2?. 



X L(y) Z Js 

e 



In order to obtain a non-trivial bound we choose 



5.17 


<£ = 111 $-, 


3 1 o 


and hence 




5.18 






* z^Tl l Sf H*) Z J ^ 

< <> s> 



y l }dv 

i Ifc f f Urf'Js . 



So that the last equation may be employed to bound 5*10 
we rewrite the previous equation as follows 



5.19 



ifVs x'^Vi/ 

8 ' z ,y D *- 



vJ!±- 





and then choose £ ^ so that the coefficients on the right 
are equal. We solve the equation 



5.20 




+ 2 Vi 



m l M 2 = 0 



and obtain 
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5.21 E 1 = 4" 1 [-2M 1 + + 8m 1 M 2 

Since 6 is to be a strictly positive number, the positive 
square root is chosen and thus 

5.22 £ 1 = \ (~M X + ^ + 2 rai M 2 ) , 



To shorten the notation we let 



5.23 



M, = TT 
-> m. Mo L 



\ (-Mj. + \[vi^ + 2m 1 M 2 ) + M x 
and bound 5-19 with 5*l8 as follows 

5-2lt | L 2 «/j 

Since a non-trival bound of the left side of the last 
equation is desired we choose the arbitrary number t ^ 
in an appropriate manner: 

Hr q 



5.25 



= 

2 M- 



Substitution of 5*25 into 5*24 yields 

5.26 iV<ij iA * S [ f l clv = IS 

which gives an estimate for the unknown boundary term of 
5.5 when this term is first bounded by 5*7* 

To bound the unknown volume term of 5*5 the following 
cases are considered. For N = 2 or N = 3 the results of 
pages 5&2 and 5&3 of [21 ] are applicable since from 5*18 
and 5.26 we have 



12 ? 



5.27 fit ftj* < ttfss i 

For N >3 a function ^ is defined by (4»l8) of [2lJ so that 



5.28 



< M 



x € D 



and 



5.29 



BC^) < o 



x 6 D 



From the divergence theorem, it can be shown that 



5-30 



tSS{&P<f8(f) - <P l B&OJJv 



« 2 SJ + Z [filin'* Vjs}^ 



+ V V 13 = } 



where 

s.31 \ = m B X ("iv" + bln ^l • 

Hence 

5-32 {- 



The desired result is now obtained by applying Schwarz's 
inequality to the unknown volume integral which we are 
seeking to bound. We obtain 
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5.33 I XT ^ 5 (q, ) ,/ V N {-XTs r V v} ^ 

and note that the first factor on the right is bounded by <, 5-32 
while the second is computable due to (4*1) of [2l] . 

Prom the above it is clear that 

5 - A ln f )-x(p)l. 



€ 






where V is given by £. 26 , the second term on the right is 
bounded by 5*33 and 

S.3S Z(p)-- Js + SSPpFjv 1 - lf(p) . 

2. A Generalized Dirichlet Problem 

In this section a pointwise bound for the solution 
of a Dirichlet problem shall be obtained where the given 
differential operator is such that the results of |2l] are 
not applicable. 

Since the problem and the techniques which are to be 
employed are very similar to those which were used for the 
Degenerate Dirichlet Problem, the notation of J2l] shall be 



129 



altered so that it more closely follows that which was 
previously developed in this paper. 

We seek a pointwise bound for the solution of 

^J( w) = (a^w .) .-qw = f(x) x6 D 

» - 1 - » J 



5.36 



w(x) 



= h(x) 



x e b 



where D and B are as defined in the first section of this 
chapter. The solution function w is assumed to be continu- 
ously differentiable in D and to have a continuous second 
derivative in the interior of a finite number of regions the 
sum of which is D. The symmetric matrix a^ is piecewise 
continuously differentiable and satisfies 



£.37 



O ta t f « 4 V 4 U.T!? 



for all real numbers (J where equality holds 

on the left if and only if all the numbers are zero. The 
bounded integrable function q = q(x) is very similar to the 
analogous function of the Degenerate Dirichlet Problem! that 
is, q has the form 

A 

5.38 q = £ I q T = 51 

i r=» I 

where 

5.39 

1=1 



and for N 3 3 
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5 . 4 ° 




Conditions similar to 4»8 a and 4*® t> may easily be obtained 
for N = 1 and N = 2. To shorten the argument, we consider 
in complete detail only the case N>3* As for the degenerate 
problem, we define m as 



So that the divergence theorem and Schwarz's inequality are 
applicable we assume that f is integrable and square integrable 
over D. 

To obtain the desired bound for the solution at an 
interior point p € D we choose an arbitrary function if = tf(x) 
which is such that J( If ) and If (x) approximate f(x) and h(x) 
as appropriate. We further require that L f be twice piecewise 
continuously differentiable in D. We introduce the function 





5.42 



V (x) = w (x) - if (x) 



and compute 





From (6.8) of [2l] we have 



5.44 



w (p)- dUp) 
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where 

s.1,5 7(r p ) - (.iJp pfl ) tJ -qP p 

with n p given by (4*2) of the work by Payne and Weinberger. 
The unknown terms on the right of 5*44 are the second and 
fourth. The final term is the first to be estimated. With 
B replacing B(y 0 ) and D replacing D(y 0 ) we note that 
equations 4*78 to 4»8l are valid for the present problem. 
Next we use 



5.46 


(I * 7/ W t j~k SScl*% x % h 

0 B 4 N q ° 0 4 * 


to write 




5.47 


I mjts « 5c ♦ ^(c-, * *■$ ♦ )su\ i% j v 




+£„' i^hyjs * # Ms 


where 




5-48 


_ -Ik 

m, = min n f n^ 

J B K 


5.49 


R 2 - min max l (xW) 2 
xeDzfB ' 


S.So 


A/ 2 

M-. = max H (f 1 ) 

D ' 


5.51 


M i}i = max { q 2 } 
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and, as indicated above, 7C is given by 4*80. To determine 
a bound for the unknown term on the right of 5*4 7, the 
divergence theorem is used as follows 






“ If %.&' *>, + ft | vh “JT<? 

0 D * 8 B 0 v 

A 

< -II + § H 5$ l5~Y^C T ^rJ^)§^-4K L ^^ Z M 

D A c b r* i + b 1 1 A 




where we have use of an appropriate form of 4*4^» Trans- 
position of £*£2 yields 

5.53 Jv 

A 

* "21 i /m r few) | (x - ^ ) ^ r ' 2 /n . /iVs t If p\*| V 

^ ^ H%k + ^ 

for any positive numbers and Since we are seeking 

a bound for the Integral of the square of the conormal 
derivative over B we are led to choose the arbitrary 
constants in the last Inequality as follows 

N^a 

5-54 «! = d-5) jjp* 
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and 

s-ss 



■a - irv 



2Mi 



= s.d-5) (o n ♦ ^ + ftMI.)- 1 

1 a o 

o 



Substitution of these results into S»i\l yields 



5-56 



$ ( §v) js < a 

6 



IS 



where Id is an appropriate computable constant which is such 

tnat 



5.57 



13 = 0 



if 



5.58 p = h = o 

For N = 2 or K = 3 we use the results of pages £62 
and 563 of [21] to estimate the unknown volume term in 5.44° 
These results are directly applicable since 5*53 provides a 
bound for 

5.59 

0 1 ^ 

If N > 3 then ~F is defined by (4»l8) of {21J so that 

5.60 ~T * M ) J *&) <0 x6D 

where 

5.61 J*(90 = (>’Jf ) ,-s^t’ . 

9 J > ± 










- 
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The divergence theorem may be used to establish that 

5.62 j(f)- ¥>*• JiV 

mvjj 1 1 
-$H*34s 

3 



and, since the right side of this equation is computable, 
it is evident that the desired pointwise bound is available 
using the techniques and the results of the paper by Payne 
and Weinberger. In the above equation is given by 5 >» 5>6 
while jy is the bound for the integral of a^J° f ^ . 

9 1 9 J 

which is obtained from 5*53 after application of 5«56>° 



CHAPTER VI 



EXAMPLE 



1. Explicit Non-Linear Mixed Problem 

In this section a bound shall be obtained for the 
solution of a non-linear boundary value problem which is of 
the type considered in the first section of Chapter III. 

We suppose that V is a three dimensional domain 
bounded by the hyperplanes y = o and y = 1 and by the surface 
whose equation is 



The boundary value problem is explicitly stated as follows: 



6.1 





y w a w aw , \-i . ? 

* tl? 1 + I?' -w- V3 : 



^W 



6.2 



3 - < 00 



N VyW 

vy(^o) = o Cx^eDtp) j |v = / 



where 



6.3 



A ^ o ® 



Prom the above we obtain 
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/Yr\ t - I 




6.4 



< 






a, = / 



v 



By inspection one observes that the differential 



operator of 6.2 is slightly different from that of 3»1« 

The remarks which are found at the conclusion of Chapter 
III suggest the reason for this difference. The additional 
tern is combined with the differential operator in order to 
indicate the improvement which follows from such an 
association. If this were not done, then 



and by close inspection of the results which are obtained 
below, it is evident that the sharpness of the error bound 
is materially diminished. 



6.4 a 



1 

Mt = - 0 + 1 

1 A 2 



Following the techniques of the earlier chapters we 



let 



6.5 




and then calculate 
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f 



p = (A + - ir) 



i ^2% “1 



6.6 





Even though CP satisfies the boundary conditions on S and 
for large A approximates the solution of the differential 
equation, it is evident that does not closely coincide 
with w throughout V since on D(o) the approximation is very 
poor. 



in the remarks accompanying equation 3*23 may be defined as 
follows: 



The set of auxiliary functions which were introduced 



6.7 



i = 1 and 2 



Thus we have 



6.8 





In what follows the solution of 6.2 is bounded at 
the point p = (0,1). Prom 3*1& we have 



6.9 
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where 










6.10 


*p = 


[4 it U-y)] _1 


exp [ 


-r 2 \4 (l-y^> -1 ] 


and 










6.11 


?( V 


^ + 


^x z 


- y 

p h y 



Since IP is such that L is zero on S s we find it 
advantageous to choose the arbitrary constants of 3*35 in 
a slightly different manner than in equations 3*38 to 3»42* 
Thus by setting 

6.12 cC 2 = oO 

and 

6.13 “ o, 

equation 3*35 may be rewritten as 

6.14 /// % f (« *, - fr - $.) SSS ? V V 

' Tf (**) jj A? As 

+ ^ 3 III O + . 



Substitution of this into 3*24 yields 



139 



6.15 




Due to the form of the above we set 



6.l6 




and then seek to determine the positive constants cC ^ and 



such a minimum problem is extremely difficulty hence in 
practice we ordinarily choose the arbitrary constants in 
some admissible manner and then seek to improve our result 
by a better choice of the approximation function (f> . For 
example we assume that 



so that the right side of 6.15 is minimized. To solve 



6.17 




and set 



6.18 




and 





From the above and 6.l6 it follows that 



6.20 



K = 0. 



Equation 6.15 may now be rewritten as 



^•21 fJ ■cl 5 ' 3 + ft czrTF) ^ - 'fS . 

We use 6.10 to obtain 

6.22 {tt^rfTds}^- {.ISlY 1 - 

and hence 
6 - 23 

vS 3 

< .3 8V 13^ 



where 13 is given by 6.21. 

To bound the first term on the right of 6.9 we 
introduce an auxiliary function which is such that the 
conclusions of Theorem 2.2 are valid. The auxiliary 
function used in this example is not identical to that 
defined by 2.39» While the original function could be 
used for the present problem, advantage is made of the 
simplicity of the differential operator to obtain a function 
which lends itself to the integrations which are required 
for the determination of the desired bound. Let 

6.24 ^ - (1-y)^ exp [-r^U-y)^ 1 ] , 



Clearly 
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6.2S for (x,y) € V 

and since 



6.26 


- 

TC^O - sx' 1 + ^ 




= ~"k 0 -q) ^Mf> * 0 



for (x,y) € V, it follows that the conclusions of Theorem 
2.2 are valid as was indicated above. For the present 
differential operator equation 3*5>7 becomes 



6.27 


MS [z<¥ (?) - 

/ 3 " > ' 'ft) 

V' 

+ ^ + 2.(«+')TP*}^ 

= ~If T H ^ ds . 


From 3*^2 


and the above we have 


6.28 


K = K ? ^ (»1 + i p ) <0 
^ A' 1 


where the 


strict inequality holds on the right due to 6.17* 



In what follows we set 



/ 



6.29 



k 2 = 0 



142 



to simplify the numerical calculations. (If greater 
accuracy were desired, one would of course take equal 



to A" 2 - 1 ) . 


Prom equations 3*57 to 3*^5 we obtain 


6.30 


J 11 If* * f - c=r)l 




^ 'fS J", j f . z t'y = ^ 



To use 3.69 it is also necessary to compute 



6.31 


/j/tJcnrCF [iKiO-d] J v 

^ ^ JIT- 0 ~ 1 R ( i-^j j J,V 




- 2 .zi R 1 



and 



6.32 


mi- i«r j*c tffjv 

V r 

= ( wT z Iff 1 0 -# * a 1 (v ( 1 ■ - 3))" '] 4 ^ 




- 0.0 R / 8 



By combining the foregoing and equation 3.69, it; follows 



that 
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6.33 



I fff V J*(0 M 

^ ■* i v 6a) 

V 



^ 0, z oy 






*.2/7 t-/.6Y 



? 



* rr ] -2]^\ 



To bound the second term on the right of 6.9s we 
could use the techniques associated with equations 3*69 
to 3*7i?> however, since 

6.34 j*( tfp) = - r p , 

we apply Schwarz’s inequality to the expression under con- 
sideration and write 

6-35 |iUu fff V« {-ft*, ^ J w\ \7w) 4 )} Jv I 

v w * 

^ 7) fff I ^ I dv 

< k { ff/tm]"f 

* 0 . 20 * a ' 1 {/?/*■' +{/,<*[ 3 + 

+ 3 . 2 V /T* 1 -^/ 7 } / . 
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The bound for the solution is completed by calculating 
the function X (p) which is given by 3. 17. One easily 
obtains 

6.36 ^(cy) +■ XT els - O 

<5 

6.37 If Cr t„ J.S - -0.06 

m P 

and 

6.38 -.IZIZR~ Z < - fSS * -22 

Prom 6.9 and the above we write as the desired bound for 
the solution 



6-39 -0.3Sv{? f -0. 2 01 + (/ f 2 A 

4- {+.2/7 f !.(,¥[ f 4* + 



-/ 



+ 0 ,o& + ,21/Z (/) 4 i)~ l 



* w(0,/) 

*+o.3M{¥ 4- 0.7b^(HR‘ l )(l.8fi l 



4 - 0 , 0(9 i-.ZZlt P 
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For example let U3 suppose that 

6.40 A - 4 

then one obtains from 6.42 the following result 

6.41 -0.95*w(0,l) s +1.09 

When bounds for the solution of more complicated 

problems are desired, the techniques are exactly the same 
as used above j however the integrals which result may be 
quite difficult to evaluate. 
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